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ABSTRACT 



For every orthonormal matrix X there exists a skew- symmetric 
matrix N such that X = (I"N)(I+N) \ provided that (X+N ) 1 exists. 

A matrix K, similar to N, can be defined for biorthonormal matrices 
U and V such that U = (I~K) (I+K) 1 , provided that (X+K) 1 exists. 
Numerical methods are presented for examination of the properties 
of K. The particular property anticipated for K, that it exhibit 
n(n-l) basic parameters inherent in biorthonormal matrices, is not 
apparent in the numerical examples derived. 
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1 . INTRODUCTION 



If A is a real symmetric matrix of order ti, it possesses n real ortho- 
gonal eigenvectors which may be normalized, as described in Appendix 1. 
The square modal matrix X‘ f of these eigenvectors, 

X' = (x x x 2 . . . x n ) (1) 

is orthonormal, that is: 

'x'X' = I = X' *X' (2) 

For the purpose of this analysis the eigenvalues of A are assumed to be 

discrete and the eigenvectors unique. 

2 

The n elements of such an orthonormal matrix are subject to the 
conditions of equation (2), of which only n(n+l)/2 are independent, leaving 
n(n-l)/2 independent parameters by which the orthonormal matrix can be ex- 
pressed. Consequently, as has been shown by Heading [3] for example, X 1 
can be written in general as 

X' = (I-N) (I+N)” 1 , (3) 

provided that (I+N) ^ exists, where N is skew- symmetric and contains the 
required n(n-l)/2 parameters. Hence any real skew-symme trie matrix de- 
fines an orthonormal matrix, provided that (I+N) is non-singular. Equa- 
tion (3) may be solved for NCX^): 

N = (I+X')“ 1 (1>X*) (4) 

The proof of equation (3) hinges upon the fact that I+N and I-N 
commute, i.e. 

'X'X' = (I-N) ~ 1 (I+N) (I~N) (I+N) 1 
- (I-N)' 1 (I-N) (I+N) (I+N)" 1 
= I 

2 

since (I+N) (I-N) = I-N 35 (I-N) (I+N) ; furthermore N must be skew- symmetric 
in order that I-T = I+N as used above. 

The discussion of this property of orthonormality in the context of 
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the eigenvalue problem is of particular significance in the reversed prob- 
lem suggested by Bell [*4]. The solution of his problem for A symmetric is 
made possible by the existence of N, as in equation (3), containing just 
the right count of parameters required for the solution. As pointed out 
in T4], the n(n+l)/2 independent elements of a symmetric A should yield 
the same total number of independent parameters in the eigenvectors and 
eigenvalues of A. This exact count does appear in the n eigenvalues plus 
the n(n-l)/2 elements of N. 

This description of the property N(X'), or N(A) for A symmetric as 
the property is restated in the preceding sentence, is background for the 
work to be reported herein. In T41 Bell has shown that an attempt to ex- 
tend his analysis to a class of unsymmetric matrix A, possessing a real 
eigenvector matrix U* , would involve the property K(U') analogous to N(X'), 
where the matrix K should express n(n-l) independent parameters inherent 
in U* , or twice as many as are inherent in X* . The development of this 
concept is presented in the next section, where the multiplicity possible 
for K is pointed out, as is the need for a definition, if possible, of 
that canonical form, , which displays the required property. 

Numerical methods are set up in the present investigation for de- 
riving K(A) , with A unsymmetric, and for exploring its properties in 
association with those of the modal matrix U 1 , and of the modal matrix 
V* which is adjoint to U* . These methods are programmed for arbitrary 
matrices of order up to n = 25. 



2 PROBLEM STATEMENT 

If A is real but unsymmetric, then its eigenvectors u^ are not 
orthogonal; but the adjoint set of vectors v^> belonging to A, satisfies 
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a biorthogonality with the u^, as described in Appendix 2 0 In fact the 
square matrices of these adjoint vectors, 



u ' u 2 . . . u n ) 

v ' "(v 1 v 2 . . . v n ) 



(5) 



satisfy the conditions of biorthonormality t 

V 1 U 1 = I = U V , (6) 

if the and are mutually normalizedo 

It should be noted here that, as before, the eigenvalues of A are 

assumed to be discrete and the eigenvectors unique; A is also restricted 

to that class for which the eigenvectors are real. 

2 

The 2n elements of such real adjoint matrices in (5) are subject 
2 

to the n conditions of equation (6), all of which are independent, indi- 
cating that both U* and V 1 can be expressed with no more than n^ basic 
parameters. As shown in [4], U 1 and V* can in fact be written as 



U' = (I-K)(I+K) 
V 1 = (I+K)(I~K) 



provided (I+K) ^ exists, where K is an arbitrary real matrix. Henc^ , by 
equations (7) any real matrix defines a pair of real adjoint or biortho- 
normal matrices U* and V 1 under the provisions statedo Equations (7) 
may be solved for K* 

K = (I+U’)" 1 (I“U') 

^ <®> 

K = (V'+I) (V'-I) 

It is readily verified that equations (7) and (8), for A symmetric, 
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are identical with equations (3) and (4), since for A symmetric: 



U' = V’ = X' 



A/ 

K = N = -K 



( 9 ) 



The proof of equations (7) follows from the commutability of (I+K) and 
(I-K) , as for equation (3). The consistency of equations (8) is demon- 
strated for the general case in f4j. 

The diagonal elements of I in equation (6) are the normalized inner 

til 

product, or dot product, of the i columns, or vectors, of U 8 and V*. 

Unlike the case for the vectors in equation (2) , this is not a unique 

specification for the dot product of u^ into itself, or of v^ into itself, 

i.e., for the "magnitudes" of u^ and v^ c For each such pair of adjoint 

vectors u^ and v^ (belonging to the same eigenvalue of A) there is a free 

condition to be chosen arbitrarily. There are n such pairs and hence n 

such arbitrary conditions 0 This suggests that the property K(U f ) can in- 

2 

volve not as many as n basic parameters, bu,t rather as few as n(n~l)o 

An additional argument given in £4l, on the basis of the eigenvalue problem 

itself, is similar to the one given in the Introduction for symmetric A, 

2 

namely: an unsymmetric matrix A has n independent elements and should 

yield just as many total basic parameters in its eigenvalues, which are n 

2 

in number, and in its eigenvectors, which must therefore have only n -n. 

As shown in Appendix 2, there are an infinite number of pairs of 
biorthonormal matrices, satisfying equation (6), belonging to a given 
matrix A. This multiplicity corresponds to the fact that the eigenvectors 
can be multiplied by any scalar and still satisfy the eigenproblem state- 
ment, coupled with the free condition on the dot product mentioned in the 
preceding paragraph. In view of equations (8), there must be a correspond- 
ing multiplicity in K(U*) for any A, or in other words in K(A) o 
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If K can be found to involve as few as n(n°l) independent parameters, 

the skew*- symmetric form which K takes with A symmetric suggests that in 

the more general case these n(n°l) parameters would be off* diagonal elements, 

and that the main diagonal would remain empty. This form of K, if it should 

occur only for one particular pair of U* and V 1 out of the infinite multi- 

plicity possible, might be designated its canonical form, K c , corresponding 

to the canonical pair, U 1 and V' . The conditions on U 1 and V 1 for the 

c c 

canonical form, if it exists, have not yielded easily to analysis for n ^ 4, 
but might become apparent from a study of numerical examples Q 

The problem which was stated for this investigation may now be outlined 
as follows: 

1. Set up the computer solution for K from a given unsymmetric A, 
using equations (8). 

This has involved the use of known techniques for diagonalizing 
A; Laguerre*s method r 5 1 was adapted. It was also necessary to de- 
velop subroutines for numerical solutions of the adjoint vectors U 1 
V 1 , when the eigenvalues by Laguerre's method were known. Care must 
be taken that the given A possess real eigenvectors. The solutions 
for U 1 and V 1 are verified for biorthonormality in equation (6). 

2. Set up the computer solution for U 1 and V 1 from arbitrary K, 
using equations (7). 

This is provided as a final subroutine to Step 1, for an over- 
all check of the solution; it is also provided as a separate program 
for the investigation of U* and V 1 from arbitrary Ko The latter pro- 
gram also includes solutions for V'V 1 , IT 1 U 1 , V 1 v and^f , U t (see 

Appendix 6) to provide for possible interpretations of the numerical 
relationships among the adjoint vectors themselves. 

The objectives of Step 1 were first to provide an overall check of 

the methods developed in T41 for A unsymmetric, and second to exhibit 

numerically some examples of the property K(A) , with particular regard for 

evidence of its basic parameters. The objectives of Step 2 were to explore 

whether evidence of the conditions on these parameters might be exhibited 
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readily in the numerical characteristics of U 1 and V' , for example in or 
among the "magnitudes’ 1 of the eigenvectors, as displayed in the elements 
of the diagonals of^U* and Tv 1 « The computer programs for these in- 
vestigations are written in Fortran 63 and can accommodate matrix orders 
of 25 or lessc 

The matrix properties which are being sought out in this investiga- 
tion are known to have hypergeometric interpretations for n greater than 
three T4]. For n less than three they represent the familiar interpreta- 
tions of Euclidean three-space. Attempts to extrapolate the latter into 
higher order hyperspace can be misleading. On the other hand, it was con“ 
sidered desirable to keep n small at the outset in order to limit numbers 
of numerical data to be examined. For this reason the numerical examples 
reported herein are for third and fourth order matrices only, 

A third computer program was formulated in order to further investig- 
ate the characteristics of the matrix K. In this case the K matrix, how- 
ever, was derived from an orthonormal matrix after the latter had been 
premultiplied by a diagonal matrix. The goal here was to observe the ef- 
fects of premultiplication in the changes on the matrix K Appendix 3 
describes in detail the background for this program. 



3o COMPUTER PROGRAMS 

Two Fortran 63 computer programs were formulated in order to satisfy 
the requirements set forth in the problem Statement above, Steps 1 and 2. 
All programs were arranged to operate with matrices of order n ^ 25. 

Appendix 5 describes the computer program named REVEIG, designed for 
Step 1, which basically is a general eigenvalue program for any real square 
matrix A where, as described in Appendix 2, 
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AU* = U'D 
7v' = V'D 

and V'U 1 =1 

The derived matrices U* and V* are then applied in equations (8) to derive 
K and It , from U fi and V 8 respectively. 

The available EIG3 subroutine used to derive the eigenvalues for REVEIG 
can be used to find real and imaginary roots of a matrix and store them for 
further use f5]. A basic assumption for REVEIG was that all matrices in- 
vestigated would be required to have real roots only* Subroutine EIG3, 
however, did not derive the eigenvectors along with the eigenvalues. A 
program for the eigenvectors was derived as an additional step, as described 
in Appendix 5. It should be noted here that due to the arbitrariness of the 
solution used for the eigenvectors, a particular restricted class of multip- 
licity is built into the program, to the exclusion of other classes which 
exist, and which might bear further investigation. Some typical test re- 
sults from REVEIG are contained in Appendix 8. 

Appendix 6 describes the computer program called REVEIGl, designed 
for Step 2, which was based on equations (7) in the Introduction. Inputs 
were made to REVEIGl from the K and If matrices derived in REVEIG, in 
order to provide a check of REVEIGl by recomputing the U* and V 1 matrices 
from which K was derived. Thereafter, other arbitrary matrices could be 
used as K inputs to REVEIGl - As a part of this program the products V 1 V* , 
V*V* , U'U 1 , U 1 U • , and V 1 U 1 are also derived. Some typical results are 
shown in Appendix 8. 

Appendix 7 describes the computer program called REVEIG2. This 
program was formulated late in the investigation in order to further explore 
the effects of matrix algebra on an orthonormal matrix. As described in 
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Appendix 7, the orthonormal matrix U 1 derived in REVEIG was premultiplied 
by a diagonal matrix G of significance in a particular case. The result 
was used as an input to derive the matrix K using equation (8) in the 
Introduction. The ultimate goal was to investigate the effects of pre- 
multiplying an orthonormal matrix in this special case, and to determine 
what characteristics, if any, exist in the resulting K matrix. The signifi- 
cance of this case is discussed in Appendix 3. 

The three computer programs are designed for an input format of 
(4E20.10). The outputs are all printed out in a (7E17.10) format. In order 
that the computer results in Appendix 8 be more easily understood, the fol- 
lowing relationships are defined: 

Computer Theory 



UP 


U' 


VP 


V' 


KT 


r 


UPT 


U 1 


VPT 


T’ 



Comment cards are included in the program decks for all the programs, that 
basically describe the individual steps involved. The program decks can be 
acquired from the Department of Aeronautics of the U. S» Naval Postgraduate 
School. 

4. RESULTS AND CONCLUSIONS 

A few of the several numerical examples run from the programs described 
are included in Appendix 8, showing typical results. It was a primary ob<- 
jective of this investigation to provide an overall numerical check of the 
methods developed in T41 for A unsymmetric , namely to derive a property 
K(A) as described in Section 2. In this objective the investigation is 
successful, and the consistency of equations (7) and of equations (8) is 
verified. 
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A second objective was to exhibit typical examples of the matrix K 
with the view that typical patterns might be evident in its properties. 

In particular the so-called canonical form dependent upon n(n~X) basic 
parameters was of interest. This form was not observed, at least not in 
the form with empty main diagonal, and was not otherwise identified. In 
this connection it should be noted that a single restricted class of multi- 
plicity has been built into REVEIG, as discussed in Section 3. 

Inputs were made into REVEIGl consisting of arbitrarily assumed 
"canonical" forms, K c , possessing zero main diagonal. It was anticipated 
that the final printout of this program, consisting of the several products 
listed in equation (2-8) of Appendix 2, might show evidence of associated 
canonical properties in U fl and V* . These properties were not observed. 

Finally a calculation with REVEIG2 was made to explore the form of K 
which results from a special case of unsymmetric matrix, as described in 
Appendix 3. No significance was found in the resulting example. 
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NOMENCLATURE 



Most terras within this thesis are defined when used; however, the 
following list of the general terras used is supplied for ready reference* 



A 



a u 

A* 

B 

D 

E 

F 

G 

I 

K 

K 

c 

n 

N 




U 

u 

c 



V 

V 

i 

w 



X 

X* 



real square matrix 
elements of A 

B ^AB ^ = a derived matrix in Appendix 3 

diagonal matrix of non-vanishing , real diagonal elements 

diagonal matrix of eigenvalues assumed to be distinct 

diagonal matrix 

diagonal matrix 

diagonal matrix 

identity matrix 

matrix property belonging to A unsymmetric 
hypothetical canonical form of K 
integer 

skew- symmetric matrix 

eigenvector (column matrix) 

modal square matrix of eigenvectors u^ 

hypothetical canonical form of U 

eigenvector (column matrix) 

modal square matrix of eigenvectors v^ 

hypothetical canonical form of V 

diagonal matrix 

eigenvector (column matrix) 

orthogonal modal matrix of eigenvectors x^ 

pseudomodal matrix (Appendix 3) 
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X ^ 55 real, discrete eigenvalues of A 

/w 

( ) = transpose of ( ) 

( ) 1 » normalized form of matrix ( ) 
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APPENDIX 1 



SIMPLE EIGENVALUE PROBLEM WITH A SYMMETRIC 

The simplest possible eigenvalue problem is defined (*1,21 by 

Axi = \ ± x ± (1~1) 

where (see Appendix 4 for definition of terms): 

A = n x n symmetric matrix 

= eigenvalue for which a solution x^ exists; there are n 
solutions assumed to be discrete 

x-£ = column matrix (eigenvector) 

The n solutions may be combined into a single statement 

AX ^ XD (1-2) 

where 

X = modal matrix of eigenvectors 
D = spectral matrix (diagonal) of eigenvalues 
It is to be noted here that post-multiplication of the modal matrix 
by a diagonal matrix does not alter the eigenpro£>lem (1-2). 

After premultiplying equation (1~2) by and also after trans~ 
posing equation (1~2) and pos tmultiplying by X: 

/V A/ 

XAX = XXD 
= DXX 

Putting XX - W and subtracting: WD - DW = 0 , i.e c W^j = 0, ij 1 j , 

since the Xi are assumed discrete 0 Therefore 

XX » W (diagonal) 

i.e., X is orthogonalo Letting X 1 = XW 2 , i.e OJ after normalization 

r* 

X’X 1 = I 

Note also that X*X* = I 

In this case X* is said to be orthonormal. 
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APPENDIX 2 



EIGENVALUE PROBLEM WITH A UNSYMMETRIC: GENERAL CASE 
If A is real but unsymmetric its eigenvectors are not orthogonal 

(self-adjoint) as in Appendix 1. However A and S' have the same eigenvalues 

(assumed to be discrete) but different eigenvectors, v^ and u^ respectively 
T2j. The n solutions for each may be combined as follows: 

AU » UD (2-1) 

AV = VD (2-2) 

where (see Appendix 4 for definition of terms): 

A =* n x n real unsymmetric matrix 

U = Modal n x n matrix of eigenvectors u^ associated with A 

V = Modal n x n matrix of eigenvectors v± associated with A 

D = Spectral (diagonal) matrix of eigenvalues 

Again as in Appendix 1, after premultiplying (2-1) by *\T, and also 
after transposing (2-2) and post-multiplying by U: 

VAU ="VUD (2-3) 

and VAU = DVU (2=4) 

rs/ 

Putting VU = W and subtracting: 

0 » WD - DW 

i.e. Wij = 0, i f j, since the are assumed discrete. Therefore 

*VU = W (diagonal) (2-5) 

or U and V are said to be biorthogonal . 

Since the basic eigenproblem 

Aui = Xif 1 ! 

is unchanged by multiplication with a scalar, it follows that any scalar 
multiple of u^ is a solution. Hence if U and V are solutions to (2-1) 
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and (2-2), respectively, so also are UE and VF, where E and F are arbitrary 
real diagonal matrices (but note that only post-mul tiplication by E and F 
is permissible). It is clear that by putting 

E “ F = W ^ (2=6) 

then U'= UE - UW - ^ 

‘I 

V 1 = VF = VW ^ 

will have been mutually normalized, i.e. 

TV = I (2-7) 

In this case U 1 and V 1 are said to be biorthonormal; they are also said 
to be adjoint. 

It will be noted that if {2-1) is satisfied, then by transposition 
and commutation it is also true that 

U'V* * UV = V'U 1 * I + V'\J' (2-8) 

It is also clear that (2-6) is not a unique formula for mutual 
normalization of U and V. The values 

E = W’ 2 , F = W 

would have done as well, or any of an infinite variety of combinations, 
i.e. , there is an infinite multiplicity in the pairs of matrices which 
satisfy equations (2-8) and which are adjoint eigenvectors of A* 
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APPENDIX 3 



A PARTICULAR EIGENVALUE PROBLEM 



A special case of unsymmetric matrix is derived from the following 
eigenproblem. Let 

AU* = BU f D (3-1) 

where (see Appendix 4 for definition of terms) 

A = real symmetric matrifc 

B = real diagonal matrix of non-vanishing diagonal elements 
U 1 ® modal matrix (real) 

D * spectral matrix (real) 

This can be converted to the form of equation (2=1), Appendix 2, since 
B ^ exists: 



CU* = U’D (3-2) 

where C = B ^A is not symmetric. Therefore, as in Appendix 2, U # is 
not orthonormal, but the adjoint V 1 exists such that V'U 1 ® I, and 
the property K(C) can be defined as in the text, equations (8) 0 

On the other hand since A and B in equation (3=1) are both symmetric 
it is possible by factoring and multiplying through by B '"to rearrange 
equation (3“1) as follows F41: 

B^U - ,B^U D 

A* X’* = X 1 * D (3-3) 



where 



~ J- 

A* - symmetric = B 2 AB 2 



(3-4) 



X 1 * = orthonormal pseudomode = B 2 U' 

D = spectral matrix, unchanged 

As in Appendix 1, since A* is symmetric then X 1 * must be orthonormal 

X'*X'* - I 
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It follows that the property N(A*) exists for X 1 *, as in the text, 
equation (4). But by equation (3-4) this property N(A*) can also be 
associated with U. Since the property K(C) is also associated with U, 
then the question is suggested of how N(A*) and K(C) may be related. This 
question can be amplified as follows* Since B, or B^, or B is arbitrary 
(diagonal), then equation (3-4) may as well be written 

U' = GX 1 * (3-5) 

where G = B 2 (3-6) 

=* arbitrary diagonal matrix 

The conversion back to the context of equation (3-1), for some given G, 
may be carried out at any time* The general question can now be stated, 
namely, given X 1 *, an orthonormal matrix with skew- symmetric property 
N(X'*), pre-multiplied as in (3-5) by any G: does the property K belong- 
ing to U 1 = GX 1 * have the canonical form sought in the text, for example 
does the main diagonal of K remain empty? 
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APPENDIX 4 



Vec tor o c . 


DEFINITION OF TERMS 

.A quantity in an n^ order system represented by a column 
matrix: 

/? 


Modal matrix . . . 


w 

.A composition of all possible eigenvectors into one square 
matrix. The columns of eigenvectors are formed side by 
side. 


Transpose ... 


.The transpose of a matrix A, written A, has for its i*^ 
row the column of A and for its column the jth 

row of A, for all rows i and columns j. X is said to be 
obtained from A by interchanging rows and columns . 


Symmetric . . . . „ . * 


.If a" = A, A is said to be symmetric. A symmetric matrix 
must be square and its elements satisfy a^j = a j-p for 
all i and j. The matrix is symmetric with respect to 
its main diagonal . 


Skew-symmetric. . 


.If the square matrix A is such that = -A, it is said 
to be skew-symmetric. In a skew- symmetric matrix, 

a.. = -a . . . In particular, the elements a., of the 
ji 11 

main diagonal must be zero. 


Spectral matrix. 


.A diagonal matrix with eigenvalues (X^) as its main 
diagonal . 


Orthogonal. ..<>.. 


.If the product of a real matrix X with its transpose is 
a diagonal matrix, W, then X is said to be orthogonal t 
XX = W. Each vector (column) of X is orthogonal to all 
other vectors (columns) of X. 


Orthonormal . . . . . 


.A non-singular matrix X is said to be orthonormal if 
— /v 

X ■ X; thus an orthogonal matrix X which has been 

normalized is said to be orthonormal: XX = I. 


Biorthogonal 


U and V are biorthogonal if ^VU = W (diagonal)* 


Biorthonormal o . . 


.If = I, then V and U are bior thonormal » 
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APPENDIX 5 



PROGRAM REVEIG 

The first program was based on the analysis of Appendix 2, to solve 
for U and V from a given A: 



AU = UD (5=1) 

AV = VD (5=2) 

where A = Square matrix up to order 25x25. 

U = Modal matrix of eigenvectors with respect to A» 

V = Modal matrix of eigenvectors with respect to Ac 
D = Diagonal matrix of characteristic eigenvalues. 

The subroutine used for the eigenvalues is the standard EIG3 subroutine, 
which solves for the roots of an n t ^ degree polynomial*. The eigenvectors 
are not supplied by this standard subroutine * It was pecessary to program 
for them a solution of the homogeneous equations A - \I = 0 for each 
eigenvalue; this was facilited by arbitrarily setting one coordinate of 
each eigenvector equal to unity. The program then tests each diagonal 
element of W = VU for algebraic sign*. If W contains negative elements, 
the sign of this element and of the respective column of U is reversed, 
before U is printed out, to permit real square roots of W in accordance 
with equation (2~6) of Appendix 2, in the next step. The change of sign 
is permissible because only direction of the eigenvector is affected 
thereby. 

Biorthonormality was then effected by determining a new U and V, 
called UP and VP in the print-out, as follows: 



where E and F both equal W 



-% 



(UP) = UE (5=3) 

(VP) = VF (5-4) 

The product was then printed out using 
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! 



the new UP and VP to check for biorthonormality . 



Following equations (8) of the text the K matrix and its transpose 
were then determined from U* and V*, respectively, using the following 
relationships 



These must, of course, be consistent derivations of K; it was partially 
the purpose of this investigation to verify this fact numerically. 

The program was tested with both a 3x3 and a 4x4 symmetric A, in 
which case UP = VP as was expected, in accordance with equations (9) of 
the text. 



K = (I+UP) 1 (I-UP) 
“k = (VP+I) 1 (vp-i) 



(5-5) 



(5=6) 
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00 P * »YOST ,S/1S/2S/E/6=51 ,5 *5000. 

TM.L,E. 

PROGRAM REVEIG t- 

OD I MENS I ON A ( 25 » 25 ) *U(25*25) ,AT(25»25) » V ( 25 » 2 5 ) »RTR(25) »RTI ( 25) ,AUX 
1 (25,25) *W( 2 5,25 ) »E( 2 5 ) »F( 25) »UP( 25,25) , VP (25,25 ) ,VPT (25,25) *FK( 25, 

225) *FKT ( 25 * 25 ) ,AUX1(25,25) ,RHS(25) 

EQUIVALENCE (UP,W> * \ 

GIVEN MATRIX A, FIND THE U MATR I X ( ASSOC I ATED EIGENVECTORS) AND D ( D I AGONAL MAT 
IX OF ASSOCIATED EIGENVALUES) USING EIG3 SUBROUTINE WHICH BASICALLY USES NTH 
EGREE POLYNOMIALS TO SOLVE THE PROBLEM. EIG3 ALSO RETAINS THE REAL ROOTS AND 
IMAGINARY ROOTS FOR LATER USE IF REQUIRED. 

] READ 999, N 
999 FORMAT (15) 

IF (N) 2,2,3 

2 STOP 

3 DO 10 I = 1 » N 

READ 998 , ( A ( I , J ) , J= 1 , N ) 

10 PRINT998, ( A( I ,J) ,J=1 ,N) 

998 FORMAT (4E20.10) 

DO 15 1=1, N 
DO 15 J=1,N 

15 AUX ( I , J ) = A ( I , J ) 

CALL EIG3 (AUX,N»N»RTR,RTI ,25) 

EIG3 WILL PRINT INTERMEDIATE RESULTS AND RTR , RT I ' 

DO 11 1=1 ,N 

IF(ABSF(RTI ( I ) J-l.E-5 ) 111,11 ,12 

12 PRINT 997 

997 FORMAT ( 29H0EIG3 SAYS EIGVAL IS COMPLEX.) 

STOP 

111 RT I ( I ) =0. 

11 CONTINUE 

REARRANGE EIGENVALUES SO THA, SMALLEST APPEARS FIRST 
FORM U MATRIX SINCE EIG3 SOLVES FOR EIGENVALUES ONLY. 

ASSUME ONE COORDINATE OF EACH SET OF EIGENVECTORS EQUALS ONE. USE EQUATION A 
LAMBDA *1=0 

CALL RFRTR(N,RTR) 

DO 1 3 I = 1 , N 
DO 14 J= 1 , N 
DO 14 K = 1 ,N 
AUX ( J » K ) = A ( J » K ) 

IF (J-K) 14,130,14 * 

130 AUX(J,K)=AUX(J»K)-RTR( I ) .. ' 

14 CONTINUE 
NM 1 =N-.l 

DO 140 K = 1 » NM1 
:40 RHS(K)=-AUX(K»N) 

CALL MATALG(AUX,RHS,NM1 ,1,0, DET *25) 

16 DO 18 J=1 ,NM1 
18 U ( J , I ) =RHS ( J ) 

U ( N ♦ I ) = 1 • 

13 CONTINUE 
PRINT 995 

'95 FORMAT (9H0U MATRIX/) 

DO 150 I = 1 , N 

50 PRINT 991, I , (U( I ,J) , J=1,N) 

CALL MATRA ( N , N , A , AT , 2 5 , 2 5 ) 
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t 



DO 19 1=1, N 
DO 19 J = 1 » N 
19 AUX ( I » J ) = A T ( I , J) 
DO SAME AS BEFORE 



USING AT 



ALSO. TRANSPOSE A AND DETERMINE 
NAL 



VICE A, VVICE U 
THE V MATRIX OF 



EIGENVECTORS 
EIG3 SUBROUTINE. 



21 

994 

120 

20 



MATRIX OF EIGENVALUES USING SAME 
CALL EIG3(AUX*N,N,RTR»RTI ,25) 

DO 20 1=1, N 

I F ( ABSF ( R T I < I ) ) -1 • E-5 ) 120,20,21 
PRINT 994 

FORMAT(29HOEIG3 SAYS EIGVAL IS COMPLEX.) 

STOP 

RT I ( I ) =0. 

CONTINUE 

REARRANGE EIGENVALUES SO THAT SMALLEST APPEARS FIRST 
CALL RERTR(N,RTR) 

DO 22 I = 1 , N 
DO 2 3 J = 1 , N 
DO 23 K=1 »N 
AUX(J,K)=AT (J,K) 

IF (JrK)23>230»23 . 

AUX (J»K)=AUX(J*K)-RTR( I ) 

'.CONTINUE 
DO 231 K = 1 » NM1 
RHS(K)=-AUX(K*N) 

CALL MATALGt AUX , RHS , NM1 , 1 *0,DET,25 ) 

DO 27 J = 1 , NM1 

V ( J , I ) =RHS ( J ) • ■ 

;v ( N , I ) = 1 . 

CONTINUE 
PRINT 992 

FORMAT < 9H0V MATRIX/) 

DO 220 1=1, N 

PRINT 991,1 , <V( I ,J) , J = 1,N) 

TEST V TRANSPOSE TIMES U TO SEE IF OFF-DIAGONALS 
ONALITY. CALL RESULT W. 

CALL MATRA(N,N»V»AUX»25»25) 

CALL MAMUL (N,N»N»AUX»U,W»25»25»25) 

PRINT 990 

FORMAT ( 9H0W MATRIX/) 

DO 28 1=1, N 

PRINT 991,1, (W( I , J) ♦ J=1,N)‘ 

FORMAT ( I 5 / ( 1X7E17.10) ) 

ASSUME W IS A DIAGONAL MATRIX 
ARRANGE W SO THAT THE DIAGONALS ARE ALWAYS POSITIVE BY 
CTIVE COLUMNS OF THE U MATRIX. PRINT U ANX W MATRICES 



AND .THE «jAME DIA 



230 
23 

231 

25 

27 

22 

992 

220 



56 

990 

28 

991 



ARE ZERO WHICH ENSURES ORTH 



CHANGING 

AGAIN. 



SIGNS OF RES 





LLL 


= 1 










DO 


50 


L = 1 , 


N 






I F ( 


W ( L 


, L ) ) 


51. 


>50,50 


51 


DO 


52 


LL= 1 


, N 




52 


U ( L L * L 


) =-u 


( LL ! 


> L ) 


— 


LLL 


= 2 








50 


CONTINUE 








GO 


TO 


( 54, 


53) • 


> LLL 


53 


PRI 


NT 


995 








DO 


55 


1=1, 


N 
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55 PRINT 991 , I , (U( I , J) , J=1 ,N) 

GO TO 56 

C TAKE MINUS SQUARE ROOT OF W AND STORE IN E AND F. 

' ) 54 DO 30 1=1, N 
WW=W( 1,1) 

EF=SIGNF( SQRTFI ABSF(WW) ) ,WW) 

EF=1./EF . ' 

E ( I ) = EF 

30 F ( I )=F> r 

C POSTMULTIPLY U BY E AND V BY F AND PRINT U*E AS UP AND V*F AS .VP. 

GALL MAMUL1 (N,N,U,E,UP ,25 ,25 ) « 

PRINT 980 

980 FORMAT ( 1 5H0U-PR I ME MATRIX/) 

DO 31 1=1, N 

31 PRINT 991,1 »(UP( I , J ) , J = 1 , N ) 

CALL MAMUL1 (N,N,V,F, VP, 25,25 ) 

PRINT 981 

981 FORMAT ( 1 5HOV-PR I ME MATRIX/) 

DO 32 1=1, N 

32 PRINT 99 1 , 1 , ( VP ( I , J ) ,J = 1,N) 

C TRANSPOSE VP 

CALL MATRA ( N, N, VP, VPT ,25,25 ) 

C MULTIPLY VPT TIMES UP TO CHECK FOR ORTHONORMALITY BY SEEING IF RESULT IS IDE) 
CITY MATRIX. 

CALL MAMUL ( N,N»N»VPT , UP, AUX, 25,25,25) 

PRINT 982 

982 FORMAT ( 32H0V-PR IME-TRANSPOSE TIMES U-PRIME/) 

C COMPUTE K AND KT MATRIX WITH GIVEN UP AND VP US.ING THE FOLLOWING EQUATIONS, K 
C (I+UP) INVERSE TIMES ( I —UP ) .AND KT= (VP+I) INVERSE TIMES ( VP— I ) 

DO 33 1=1, N ‘ 

33 PRINT 991,1 ,(AUX( I ,J) , J = 1 , N ) 

DO 40 1=1, N 

DO 40 J = 1 , N 

IF ( I-J) 413,414,413 

414 AT ( I , J ) =1.+UP( I ,J ) 

GO TO 40 

413 /VT ( I ,J) =+UP( I ,J) 

40 CONTINUE 

C USE AT AS TEMPORARY STORAGE 

CALL MATALGI AT , AUX 1 , N , N , 1 , DET ,25) 

DO 41 1=1, N 
DO 41 J=1 ,N 

41 AUX ( I , J ) = -UP ( I , J) 

DO 42 1 = 1, N 

42 AUX (1,1) = AUX ( I , I ) +1 • 

CALL MAMUL (N»N»N»AUX1*AUX*FK»25*25*25 ) 

PRINT 983 

983 FORMAT (9H0K MATRIX/) 

DO 43 I = 1 , N 

43 PRINT 991,I,(FK(I,J) *J = 1»N) 

DO 44 I =1 ,N 
DO 44 J = 1,N 
IF (I-J) 415,416,415 

416 AT ( I ,J)=1.+VP< I , J) 

GO TO 44 

415 AT ( I , J ) =VP( I ,J> 
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44 CONTINUE 

gSE AS TEMPORARY STORAGE 

CALL MATALG( AT,AUX1 ,N*N»1 *DET , 25 ) 

DO 45 1=1, N 

DO 45 J= 1 »N t 

45 AUX ( I ,J)=VP( 1 »J) 

DO 46 1=1, N 

46 AUX ( I , I ) = AUX ( I » I ) -1 • 

CALL MAMULI N,N,N, AUX 1- , AUX , FKT , 2 5,25,25) 

PRINT 984 

§84 FORMAT ( 10H0KT MATRIX/) 

DO 47 1=1, N 

47 PRINT 991,I,(FKT( I»J) ,J=1,N) 

GO TO 1 

END 

SUBROUTINE EIG3(A,N,M»RTR,RTI,NQ) 

474 10152 EIGENVALUES OF REAL MATRICES 
051/1 EIGENVALUES OF NON-SYMMET R I C MATRICES 
D I MENS I ON A ( NQ»NQ) , NC ( 1 00 > , RTR ( M ) , RT I ( M ) 

CALL OVFSET 
TRACE= A (1,1) 

DO 10 I = 2 , N 

10 TRACE=TRACE+A( 1,1) 

WRITE OUTPUT TAPE 6,4*TRACE 
CALL TRING(A,l»E-7 ,N,NC, NQ ) 

TRACE=A (1,1) 

DO 11 1=2, N 

11 TRACE = TRACE+A( 1,1) 

WRITE OUTPUT TAPE 6, 4, TRACE 
NU = 0 
NV = 0 

13 IF (NV-N) 14,12,14 

14 NV = NV+ 1 
NU = NV 

16 IF (NCI NV) ) 15,17,15 

15 NV=NV+1 
GO TO 16 

17 IF (NV-NU) 19,18,19 , 

18 | RTR(NU) =A(NU,NU ) 

RT I ( NU ) = 0 . 

WRITE OUTPUT TAPE 6 , 5 , RTR ( NU > , RT I ( NU ) 

GO TO 13 ‘ ;■ 

19 IF ( NV-NU-1 )20,21,20 
20 NP=XM I NOF ( M , NV ) 

CALL LAGER( A , 1 . E-4 , NP , NU , NV , RTR , RT I , NQ ) 

GO TO 13 

21 RR=.5*(A(NU,NU)+A(NV,NV) ) 

ElsRR#*2-A(NU»NU)*A(NV»NV)+A(NU»NV)*A(NV»NU) 

S=SORTF (ABSF (El ) ) 

IF ( El ) 22 ,23 ,23 
23 RTR ( NU ) =RR+S 
RTI (NU) =0. 

RTR(NV) =RR-S 
RT I ( NV ) =0 • 

25 WRITE OUTPUT TAPE 6 , 5 , RTR ( NU ) , RT I ( NU ) , R TR ( NV ) , RT I ( NV ) 
GO TO 13 25 

l 



E I G 3 0 0 ( 
E I 6 3 0 0 ( 

i 

ADDED 
E I G300( 
E I G 3 0 0 f 
E I G300( 
E I G 3 0 0 ( 



EIG3 00!- 
EIG300: 
E I G300 : 
EIG300: 
EIG300! 
E I G 3 0 D ] 
E I G 3 0 0 j 
EIG3001 
E I G 3 0 o : 
EIG300: 
EIG3002 
EIG3002 
EIG300; 
E I G 3 0 0 i; 
EIG300I 
EIG300J 
EIG300I 
EIG300; 
EIG300I 

EIG300: 
EIG3003 
EIG300- 
EIG3001 
EIG3003 
E I G 3 0 0 : 
EIG300: 
EIG3002 
EIG3002, 
EIG3002 
E I G300^ 
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RTR (NU) =RR 
RT I ( NU ) =S 
RTR (NV) =RR 
RT I (NV) =-S 
GO TO 25 
12 X=0 • 

C CALL FPOLD 

t DO 24 J=1,M 

2(f X = X + RTR ( J ) 

WRITE OUTPUT TAPE 6,6,X 
). RETURN 

4 FORMAT( 1H048X ♦ 7HTRACE =E16.8) 

5 FORMAT ( 11HOEIGEN VALUE 12X,2E20.8) 

6 FORMAT! 1H03 5X » 20HSUM OF EIGENVALUES =E16.8) 

END 

C051/2ALMOST TRIANGULAR (HESSENBERG) SUBROUTINE 
SUBROUTINE TR ING ( A , EPS , N , I NT * NO) 

DIMENSION A ( NQ » NQ ) » INT(NQ) 

.WRITE OUTPUT TAPE 6*1 

N 1 =N-1 

N2=N-2 

DO 21 J-1.N1 
S=ABSF(A( J,J+1) ) 

J1=J+1 
J2 = J + 2 
L= J 1 

N J 1 = N- J 1 

IF ( N J 1 ) 15,15,6 
6 DO 12 K=J2 » N 

T=ABSF ( A ( J , K ) ) 

IF(T-S) 12,12,11 

11 L = K 
S = T 

12 (CONTINUE 

I F ( L-J 1)13,15,13 

13 DO 131 K.= 1 * N 
T = A ( K , J + l ) 

A(K,J+1)=A(K,L) 

131 A ( K » L ) =T 

14 DO 141 K = 1 , N 
1 T=A ( J+l ,K ) 

' A ( J+l , K ) = A ( L , K ) 

141 A ( L ,K ) =T 

15 I F ( S-EPS*M I N 1 F ( ABSF ( A ( J , J ) ) , ABSF ( A ( J + l , J+l ) ) ) ) 16,16, 17 

16 L=0 

N J 1 =0 
GO TO 181 

17 T = A ( J » J + 1 ) 

DO 18 K=J2 , N 

18 |A(J,K)=A(J,K)/T ' ' . 

181 IDO 20 1=1, N 

•_ ' i M=MIN0F( J, 1-2 ) 

U=0 • 

IF ( N J 1 ) 19,19,7 
7 DO 8 K=J2,N 

8 ‘ U=U+A(K,I )*A( J,K) 



E I G3or 
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9 IF (M) 20,20,9 

9 DO 10 K=1,M 

10 U=U-A(K,I )*A< J+1,K+1 ) 

20 A( J+l, I ) = A ( J+1,I )+U 

21 I NT ( J ) =L 

22 I NT ( N ) =0 
RETURN 

FORMAT ( 1H048X » 22HALMOST TR I ANGULAR FORM)’ 

END 

0$1 /3LAGUERRE METHOD 

SUBROUTINE LAGER ( A , EPS , N1 , NU , N , RTR , RT I ,NDI ) 

DIMENSION A ( ND I ,NDI ) *P ( 6 » 1 01 ) ,RTR(NDI ) ,RTI (NDI ) ,0(6) 
WRITE OUTPUT TAPE 6,1 
0NCE=0 • 

BL1=1. 

NUQ=NU-1 
LLY = 0 
DELOLD= 1 • 

ROLD=l • 

EGSUM1 = 0 • 

EGSUM2=0. N • 

CALL FPTEST(Z) 

DO 10 L = 2 » 6 

10 P ( L , NU ) =0 . 

NU1=NU+1 

CUP=0. 

56 11 J=NU 1 , N 

11 CUP=CUP+ABSF( A( J-l , J ) ) 

CUP=CUP/FLOATF( N-NU) 

CAP=0. 

P ( 1 , NU ) = 1 . 

fIND TRACE OF H AND H SQUARED 

SPUR 1=A ( NU , NU ) 1 

SPUR2=A(NU,NU)**2 
DO 13 J =NU 1 , N 
SPUR1=SPUR1+A( J,J) 

13 SPUR2=SPUR2+A( J,J)**2+2.*A( J-l ,J)*A( J, J-l ) 

( 

INITIAL ITERATE FROM INFINITY 

14 S1R=EGSUM 1-SPUR 1 
S2R=SPUR2-EGSUM2 
F1=N-NUQ 

>IF< ABSF(S1R)+ABSF(S2R)-1.E-7*CAP) 15, .15, 16 ' 

15 XBAR=CUP 
YBAR=0 . 

GO TO 23 

16 F2 = FI- 1 • 

DR=F2* ( F1*S2R-S1R**2 ) 

L ER=SQRTF ( ABSF ( DR ) ) 

IF(DR) 17,18, 18 
17 XBAR=-2«*S1R/F1 
YBAR=2 • *ER/F1 
i GO TO 23 
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EIG301: 

EIG301 • 

EIG301: 

EIG301! 

EIG301! 

EIG301- 
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18 Y3AR=0. 

F2=S I GNF ( 1. .SIR ) 

IF (SIR) 2 2. 20.22 
20 F2=0 • 

- 22 XBAR=-(S1R+F2*ER)/F1 



EVALUATE POLYNOMIAL AND DERIVATIVES 



23 IF (ABSF(YBAR)-AB$F(X8AR)*l.E-6 124,25.2 5 

24 YBAR=0 • 

25 M=6 

I F ( YBAR)27. 26*27 

26 M=3 

27 DO 34 K=NU ♦ N 

T = — A ( K . K + 1 ) , 

DO 34 L = 1,M 

S=S I GNF ( l.,3»5-FLOATF(L) ) 

M1=L+3*XFIXF(S) 

R=-XBAR*P( L,K)+YBAR*S*P(M1 , K ) -FLOAT F< XMODF ( L-l ,3 ) ) *P ( L-l . K ) 
DO 28 J=NU ♦ K 

28 R=R+P(L,J)*A(K, J) 

CALL OVFTST (Z) 

IF(Z)29. 32.29 

29 Z = 0 • 

'P( 1 »NU)=1.E-10*P( l.NU) 

IF(P(1,NU) )30. 30.27 

30 F = FLOATF ( K-NU )/ FLOAT F ( N-NU + 1 > 

WRITE OUTPUT TAPE 6 ♦ 2 * XBAR , Y3AR , P ( 1 , NU > , F . PNC E 

XBAR=XBAR#F 

YBAR=YBAR*F 

P(1»NU>=1«0 

GO TO 27 . 

32 I F ( N-K ) 33 , 33.34 

33 T= 1 • 0 

34 P(L.K+1)=R/T 



SCALE DOWN 



39 

35 



36 



DO 39 K= 1 » 6 
B ( K ) =0 • 

DO 35 J= 1 . M 
B(J)=P( J.N+1) 

G1 = ABSF ( B ( 1 ) ) +ABSF ( B ( 4 > ) 

G2=ABSF(B(2) )+ABSF(B(5) ) 

G3 = A8SF ( B ( 3 ) ) + ABSF(B(6) ) 

WRITE OUTPUT TAPE 6 ♦ 2 . XBAR , YBAR ,G1 .G2 *G3 
FORMAT (8H I TER ATE20X . E 1 5 . 8 . 5X , E 1 5 . 8 . 8X , 3E 1 5 . 4 ) 
D=ABSF ( B ( 1 ) ) 

DO 36 K = 2 »M 
D=MAX1F(D.ABSF(B(K) ) ) 

CALL SCALE ( D * B ♦ M) 

IF (Gl) 41,41*43 

REMOVE KNOWN ROOTS 
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EIG301 

EIG301 
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EIG301 
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ADDED 
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EIG301 

EIG301 

EIG301 

EIG301 
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EIG301 

EIG301 

EIG301 
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EIG301 

EIG301 

EIG301 

EIG301 

EIG301 
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E IG301 
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EIG301 
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EIG302 
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Q1 I =0. 

Q2R = 0. 

Q2 I =0 • 

I'F(NUQ-NU) 19.21*21 
21 DO 44 J =NU * NUQ 
D1=RTR( JJ-XBAR 
D2=RTI ( JJ-YBAR 
D=Dl*-*2+D2**2 
D1=D1/D 
D2=-D2/D 
Q1R=Q1R+D1 
011=01 I+D2 
Q2R=Q2R+D1**2-D2**2 

44 Q2I =021+2. *D1*D2 

FIND SI AND S2 

19 T 1R = B ( 2 ) / 3 ( 1) 

T 1 1 =0 . 

T2R=B(3)/3(2) 

T2I=0. 

IF(YBAR)45, 46*45 

45 D1=B(1)**2+B(4)**2 
D2=B(2)**2+B(5)**2 

T1R=(B( 2 )*B( 1 )+B( 5 >*B(4) ) /D1 
T1I=(B(5)*B(1)-B(4)#B(2) ) / D 1 
T2R=(B(3)*B(2)+B(6)*B(5) )/D2 
T2I=(B(6)*B(2)-B(5)*B(3))/D2 

46 S1R=T1R+Q1R 
$1 1 =T1 1 +Q1 1 

S2R=T1R*(T1R-T2R)-T1 I*(T1I-T2I J-Q2R 
S2I=T1R*(T1I-T2I ) +T1 I* ( T1R-T2R )-Q2 I 

F.IND THE NEXT ITERATE 

LLY=LLY+1 

D=ABSF( XBAR )+ABSF-( YBAR ) 

IF(l.E+7-D*(ABSF(SlR)+ABSF(SlI ) ) ) 41,41,42 

41 MARK=1 

GO TO 100 

42 G=N-NUQ 

48 IF(YBAR-ABSF(X) >50,50,49 

49 S1I=S1I+1./(2.*YBAR) 

S2R=S2R+1./(4.*YBAR**2 ) 

G=G- 1 . 

50 IF (811)65, 65,66 

65 H=.5*(G-2.) 
pO TO 67 

66 H=G-1 . 

67 DR=H*( G*S2R-S1R**2+S1 1**2 ) 
DI=H*(G*S2I-2.*S1R*S1I ) 

IF ( DI ) 53*51 *53 

51 El =0. 

ER = SQRTF( ABSF(DR) ) 

I F ( DR ) 52 , 54 , 54 

52 El =ER 
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EIG3024 
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EIG3025 
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EIG3025 
EIG3025 
EIG3025 
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ER = 0 . 

GO TO 54 

53 CALL CXSQRT (DR*DI » ER » E I ) 

54 IF(S1R*ER+S1I*EI )55>56»56 

55 ER=-ER 
El— El 

56 D1 =S 1R+ER 
D2 = S1 I +E I 

D = D1**2'+D2 m ’*2 

X=-G*D1/D 

XBAR=XBAR+X 

Y=G*D2/D 

YBAR=YBAR+Y 

DELNEW= ASSF ( X ) +ABSF ( Y ) 

RNEW=DELNEW/DE'LOLD 
D=ABSF ( XBAR ) +A8SF ( YBAR ) 

TEST FOR LINEAR CONVERGENCE 

IF(LLY-3)62,62,57 

57 IF (DELNEW-MAX1F(3.*DEL0LD,.5*D))571, 571,570 

570 IF (BL1) 571*571*572 
572 DELOLD=CAP 

ROLD=3 • 

IF (LLY-15) 14,14,100 

571 I F ( RNEW- • 7*ROLD ) 62,58*58 
58 MARK=3 

; IF (DELNEW-.001*EPS*CAP ) 70,-59,59 

59 IF(BL 1)61,61* 60 

60 XBAR=XBAR-X 
YBAR=YBAR-Y 
BL1=0. 

GO TO 48 

61 BL 1 = 1 • 

•GO TO 63 

; j 

TEST FOR AN EIGENVALUE 

62 I F ( DELNEW-EPS*MAX1F( D, .001*CAP ) ) 64,64,63 
63 DELOLD = DELN EW 
•ROLD=RNEW 

IF (LLY-15) 23,23,100 ‘ .. ' 

DO WE HAVE A COMPLEX APPROACH TO A REAL ROOT 

64 MARK =2 

70 BL 1=1 

IF (YBAR )71 , 100,71 ' 

71 IF(G2*ABSFCYBAR)-G1)72,100,100 

72 < IF(ONCE)73,73,100 

73 X=0. 

0NCE=1 • 

! YBAR=0 • 

| GO TO 63 

- WE ACCEPT (XBAR, YBAR) AS A ROOT 



F I G 3 0 2 
FIG302 
EIG3C2 
EIG302 
EIG302 
EIG302 
• EIG302 
EIG302 
EIG302 
EIG302 
EIG302 
EIG30? 
EIG302 
E IG302 
EIG302 
EIG302 
EIG302 
EIG302 
E I G 30 2 
E IG302 
E IG302 
EIG302 
EIG302 
EIG302 
E IG302 
EIG302 
EIG302 
EIG302 
E I G 30 2 
EIG302 
EIG302 
. EIG302 
EIG302 
EIG302 
EIG302 
EIG302 
EIG303 
EIG303 
EIG3C3 
EIG303 
E I G 3 0 3 
E IG303 
EIG303 
EIG303 
EIG303 
E IG303 
E IG303 
FIG303 
FIG303 
EIG303 
EIG303 
EIG303 
EIG303 
E IG303 
EIG303 
EIG3C3 
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100 NUQ=NUQ+ 1 

RTR ( NUQ ) =XBAR 

IF(ABSF(YBAR)-.001*A8SF(XBAR )) 74,74,75 

74 YBAR=0 • 

75 IF(NUQ-NU) 9,76, 9 

9 I F ( RT I ( NUO-1 )) 76,76,77 

76 YBAR=ABSF( YBAR) 

RT I ( NUQ ) =YBAR 
GO TO 78 

77 RTI (NUQ)=-ABSF( YBAR) 

78 WRITE OUTPUT TAPE 6 , 3 , RTR ( NUQ ) , RT I ( NUQ ) ♦ LLY , MARK 

3 FORMAT! 1 1H0 E I GENVA LU E 1 2X , 2 E20 . 8 , 12X , I 3 , 1 7H I TERAT I ONS , TEST II//) 
LLY = 0 

CAP=MAX1F(D,CAP) 

DELOLD= 1 • 

ROLD- 1 . 

EGSUM1=EGSUM1+RTR ( NUQ ) 

EGSUM2 = EGSUM2 + RTR ( NUQ ) **2“RT I ( NUQ ) **2 
I F ( NUQ-N 1 ) 80 » 10 1 » 1 01 

80 IF ( YBAR ) 83 , 84 , 81 

81 YBAR=-YBAR 
GO TO 23 

83 I F ( NUQ-NU >31,84,31 

31 RTI ( NUQ— 1 ) = • 5* (RTI (NUQ-1 )-RTI (NUQ) ) 

RTI ( NUQ ) = -RT I (NUQ-1 ) 



A NEWTON ITERATE TOWARDS NEXT ROOT 

{ 

84 ONCE = 0 « 



Z = 0. 



I F ( (ABSF(Q1R)+ABSF(Q1I ))* D- 10000. )85»85, 14 
5 IF ( ABSF ( EG SUM 1-S PURI ) +ABSF ( EGSUM2-SPUR2 )-l « E-5*CAP >15,15,86 
86 DR=B(3)+2.*(B(2)*Q1R-B(5)*Q1I ) 

DI=B(6)+2.*(3(2)*Q1I+B(5)*Q1R) 

D2=DR**2+DI **2 



XBAR = XBAR-2 .*( DR*B ( 2 ) + DI.*B ( 5 ) ) / D2 
YBAR=ABSF ( Y 3AR-2. * ( DR*B ( 5 ) -D I *B ( 2 ) ) / D2 ) 

GO TO 23 
101 RETURN 

1 FORMAT ( 1H050X , 19HLAGUERRE I T ER AT I ONS /3 1 X , 9HR EAL 
1ART22X,1HP11X,7HP PR I ME6X , 1 1HP DBL PRIME) 

END 

>1 / 5 COMPLEX SQUARE ROOT 

SUBROUTINE CXSQRT ( A , B , X , Y ) 

E = MAX1F ( ABSF ( A ) ,A3SF(B) ) 

F=F*SQRTF( (A/F)**2+(B/F)**2) 

I F ( A ) 1 , 1 , 2 



PART10X, 10HIMAG. 



1 Y=SQRTF ( ( F - A ) * • 5 ) 

X=.5*B/Y 

I F ( X ) 4 , 3 , 3 
A X=-X 
Y = -Y 
GO TO 3 

2 X=SQRTF ( ( F+ A ) # . 5 ) 31 

Y= . 5*B/X 



EIG3032 
EIG3032 
E IG3032 
• EIG3032 
E I G 3 0 3 2 
EIG3032 
! EIG3032 
EIG3032 
EIG3032 
EIG3032 
EIG3033 
EIG3032 
EIG3032 
EIG3032 
EIG3032 
EIG3032 
EIG303: 
EIG3032 
E I G 30 32 
EIG3032 
EIG303A 
EIG3034 
EIG303A 
E I G 3 0 3 A 
EIG303A 
EIG303A 
EIG303A 
EIG303A 
EIG303A 
EIG303A 
E I G3035 
EIG3035 
EIG3035 
EIG3035 
EIG3035 
EIG3035 
EIG3035 
EIG3035 
E I G 30 35 
E I G303 c 
PE I G 30 3 6 
EIG3036 
E I G 30 36 
EIG3036 
EJG3036 
EIG.30 36 
E I G 30 36 
EIG3036 
EIG3036 
E I G 30 36 
EIG3037 
EIG3037 
EIG3037 
EIG3037 
EIG3037 
EIG3037 



3 RETURN 
END 

SUBROUTINE SCALE (D»B*M) 

DIMENSION B ( 6 ) 

5 = 2. ** ( X I NT F ( LOGF (D) /. 69314718056) + 1) 
DO 1 I = 1 »M 
1 B ( I ) = B( I ) /S 
■RETURN 
END 

IDENT 
ENTRY 

OVFSET SLJ 
NOP 
CALL 
11 

CALL 
ENA 
STA 
SLJ 

OVFLU SLJ 
LDA 
STA 

'• CALL 
■ . 11 
CALL 
SLJ 

OVFTST SLJ 
. LDA 
ARS 
SAL 
INA 
SAU 

■ I 1 LDA 
f STA 
ENA 
STA 
SLJ 
END 

SUBROUTINE MAT A LG ( A , X , NR » NV , I DO , DET , NACT ) 
DIMENSION A (NACT, NACT) , X ( NACT , NACT ) 

I F ( I DO ) 1 , 2 , 1 

1 DO 3 1=1, NR 
DO 4 J=1,NR • 

4 X( I ,J) =0.0 
3 ! X ( I ♦ I) = 1 . 0 

'nv=nR 

2 DET= 1 . 0 
NR 1 =NR- 1 
DO 5 K=1 »NR1 
IR1 =K+1 

i P I VOT =0.0 
DO 6 I =K , NR 
Z= ABSF ( A ( I » K ) ) 

IF(Z-PIVOT) 6,6,7 
7 P I VOT=Z 
I PR= I 



OVFSET 

OVFSET, OVFTST,OVFLU 

SELECT * 

OVFLU 

ERROR* 

0 

= SZ I G 
OVFSET 

-K-# 

= D 1 • 

ZIG 

SELECT* 

OVFLU 

ERROR* 

OVFLU 

* 

24 
* + 2 
1 

OVFTST 

ZIG 

0 

ZIG 

OVFTST 



EIG302 

EIG303 



0000 
0000 
0000 
0000 
0000 
0000 
0000 
no oo 
noon 
0000 
0001 
0001 
000 1 
0001 
0001 
0001 
0001 
0001 
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6 CONTINUE 

IF'(PIVOT) 8,9,8 
9 DET = 0 . 0 
RETURN 

8 IF(IPR-K) 10,11*10 

10 DO 12 J = K »NR 
Z = A( I PR , J ) 

A( IPR,J)=A(K,J) 

12 A ( K , J ) =Z 

DO 13 J = 1 , NV 
Z = X< I PR » J ) 

X( I PR * J ) =X ( K * J ) 

13 X ( K » J ) =Z 
DET=-DET 

11 DET = DET*A ( K * K ) 

PIVOT =1«0/A(K»K) 

DO 14 J= I R 1 * NR 
A(K,J)=A(K»J)*PIVOT 
DO 14 1 = I R 1 , NR 

14 A< I , J) =A( I ♦ J) -A( I ,K) *A (K» J ) 

DO 5 J= 1 * NV 

IF (X(K*J> ) 15,5,15 

15 X(K,J)=X(K, J)*PIVOT 
DO 16 I = I R 1 , NR 

16 X( I»J)=X( I,J)-A(I,K)*X(K»J) 

5 CONTINUE 

IF ( A ( NR »NR ) ) 17,9,17 

17 DET=DET*A(NR,NR) 

PIV0T=1.0/A(NR,NR) 

DO 18 J=1,NV 
X(NR,J) =X( NR , J)*PIVOT 
DO 18 K= 1 »NR 1 

I =NR-K 
SUM=0 • 0 
DO 19 L = I , NR 1 

19 SUM=SUM+A( I ,L+1 )*X(L+1 ,J) 

18 X( I , J)=X( I , J ) -SUM 
END 

SUBROUTINE MATRA ( M ,N ,A ,B ,MD ,ND ) 
DIMENSION A ( MD , ND ) » B ( MD »ND ) 

DO 3 1=1, M 

DO 3 J = 1,N ' •• 

3 B ( I ,J)=A( J, I ) 

! RETURN ' 

END 

SUBROUTINE MAMUL ( M ,N , L , A , B , C , MD , ND , LD ) 
DIMENSION A(MD,ND) ,B(ND,LD) ,C(MD,LD) 

DO 2 I = 1 ,M 
DO 2 J = 1 , L 



S=0 • 

DO 1 K = 1 ,N 

- 1 S=S+A ( I ,K)*B(K, J) 

2 C ( I » J ) =S 



RETURN 

END 

SUBROUTINE 



MAMUL1(M,N,A,E,B,MD,ND) 



0001P 
0001 ° 
00020 
00021 
00022 
0002 3 
0002 - 

00025 

00026 

00027 

00028 

00029 

00030 

00031 

00032 

00033 

00034 

00035 

00036 

00037 

00038 

00039 

00040 

00041 

00042 

00043 

00044 

00045 

00046 

00047 

00048 

00049 

00050 

00051 

00052 

00053 

00054 

00055 
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DIMENSION A ( MD ♦ ND ) * E ( ND ) * 3 ( MD » ND ) 
DO 1 I = 1 * M 
DO 1 J = 1 * N 
1 B( I ,J)=A( I * J)*E( J) 

RETURN 

END 

SUBROUTINE RERTR ( N *RTR ) 

DIMENSION RTR ( 25) 

NM 1 =N- 1 
DO 12 I = 1 * NM 1 
T P = I + 1 

DO 12 J= I P » N 

I F ( RTR ( I ) -RTR ( J ) ) 12*12,13 
13 TFMP=RTR ( I ) 

RTR ( I ) =RTR ( J ) 

RTR ( J ) = TEMP 
12 CONTINUE 
RETURN 
END 

END 

FINIS 

-EXECUTE * 3 « 

4 



i 

{ 

X 




7 
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APPENDIX 6 



PROGRAM REVEIGl 

The second program, REVEIGl, was formulated for a maximum 25x25 
input based upon equations (7) of the text: 



UP - (I-K) (I+K)" 1 


(6-1) 


VP = 


(6-2) 



The input is any real square matrix K. 

In addition, various operations with UP and VP are printed out for 
investigation of any obvious relationships, particularly among the diagonal 
elements; these are V'V 1 , V'V', U'U' , U'u'' , and^'U'. 

The K outputs of REVEIG were tested as inputs to REVEIGl in 3x3 and 
4x4 matrices, and in both cases the UP and VP compared with the UP and VP 
in REVEIG. V*U* always resulted in the identity matrixo When K was skew- 
symmetric, UP and VP were equal as was expected. 
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-CQOP, » YOST *S/1S/2S»3»5000. 

-FIN »l »e. 

PROGRAM REVEIG1 

’ 00 I MENS ION UP ( 25 * 25 ) *VP(25»25) »FK(25»25) *FKT(25*25) ,VPT( 25,25) , 

1UPT ( 25,25 ) ,AUX( 25,25 ) ,AUX1 (25,25)»T(25,25)‘ \ 

C GIVEN A K MATRIX FIND ITS ASSOCIATED U PRIME MATRIX USING THE FOLLOXING E0UA1 
C QN. UP=( 1-0 ( I+K) INVERSE 
>, 1 READ 899, N 
IF (N) 2,2,3 

2 STOP 

3 PRINT 870, N 

870 FORMAT (15H1K MATRIX, N= I 3 ) 

8.99 FORMAT ( 15 ) 

DO 48 1=1, N 

READ 898, (FK( I , J) ,J=1 ,N) 

48 PRINT 898, (FK( I ,J) ,J=1,N) 

898 FORMAT (4E20.10) 

C STORE ( I+K ) IN AUX 
DO 49 I = 1 ,N 
DO 49 J = 1 , N 
IF(I-J) 490,491,490 
491 AUX (1,1 ) =FK ( I , I )+l. 

GO TO 49 

490 AUX ( I , J ) =FK ( I , J ) 

49 CONTINUE 

C USE AUX AS TEMPORARY STORAGE 
C INVERSE OF (I+K) IN AUXl 

CALL MATALG(AUX,AUX1,N,N,1,DET ,25) 

C STORE ( I-K) IN T 
DO 50 1=1, N 
DO 50 J=1,N 

50 T( I , J ) =-FK ( I , J) 

DO 51 1=1, N 

51 T( I , I ) =T( I , I ) + l . 

C PREMULTIPLIES (I+K) INVERSE BY (I-K) 

CALL MAMUL (N*N,N,T»AUX1 ,UP ,25,25,25) 

PRINT 897 

897 FORMAT ( 1 5H0U-PR I ME MATRIX/) 

DO 52 1=1, N 

52 PRINT 896, I , ( UP ( I , J) ,J=1 ,N) 

896 FORMAT ( I 5 / ( 1X7E17. 10) ) 

C TRANSPOSE THE K MATRIX AND FIND ITS ASSOCIATED V MATRIX USING THE FOLLOWING E 
C UATION. VP=( I+KT) ( I-KT) INVERSE 
CALL MATRA( N,N,FK,FKT ,25,25) 

C STORE .( I+KT ) IN AUX 
DO 53 1=1, N 
DO 53 J=1,N 
AUX ( I , J ) =FKT ( -I , J ) 

C STORE ( I-KT ) IN AUXl 

AUXl ( I , J)=-FKT( I » J ) 

IF(I-J) 53,54,53 
54 AUX (1,1) = AUX ( I , I ) + 1 . 

AUXl ( I , I ) = 1 • -FKT ( 1,1) 

53 CONTINUE 

C STORE INVERSE OF (I-KT) IN T 

CALL MATALG(AUX1 ,T ,N,N»1 »DET ,25 ) 



) 



PREMULTIPLY (I-KT) INVERSE BY (I+KT) 

CALL M AMUL (N»N»N»AUX,T »VP»25»25»25) 

PRINT 895 

$95 FORMAT ( 15HOV-PR IME MATRIX/) 

DO 55 1=1, N 

55 PRINT 896, I, (VP(I*J), J=1*N) 

PRINT OUT ALSO V PRIME TRANSPOSE TIMES V PRIME* V PRIME TIMES’ V PRIME' TRANS Pf 
E. SAME FOR U PRIME AND U PRIME TRANSPOSE, 

894 FORMAT ( 9H0VPT X VP / ) 

893 FORMAT ( 9H0VP X VPT / ) 

892 FORMAT ( 9H0UPT X UP / ) 

891 FORMAT ( 9H0UP X UPT / ) 

CHECK TO SEE IF V PRIME TRANSPOSED WHICH PROVES ORTHONORMALITY. 

890 FORMAT (9HOVPT X UP/) 

CALL MATRA ( N ,N , VP , VPT ,2 5 , 25 ) 

CALL MATRA (N,N, UP, UPT, 25*25) 

CALL MAMUL ( N ♦ N , N , VPT ♦ VP , T , 2 5 ♦ 25 , 2 5 ) 

PRINT 894 
ASSIGN 56 TO I SW 
GO TO 60 

56 CALL MAMUL ( N , N , N , VP , VPT , T , 2 5 , 25 , 2 5 ) N 
PRINT 893 

ASSIGN 57 TO I SW 
GO TO 60 

57 CALL MAMUL ( N , N , N , UPT , UP * T , 2 5 , 25 , 2 5 ) 

PRINT 892 

ASSIGN 58 TO I SW 
GO TO 60 

58 CALL MAMUL ( N,N,N,UP,UPT ,T ,25*25,25 ) 

PRINT 891 

ASSIGN 59 TO I SW 
GO TO 60 

59 CALL MAMUL (N*N,N,VPT , UP , T ,25* 2*5, 25) 

PRINT 890 

ASSIGN 71 TO ISW ' . 

60 DO 61 1=1, N < 

61 PRINT 896, I , ( T ( I , J ) ,J=1»N) 

GO TO ISW 

71 GO TO 1 
END 



SUBROUTINE MATALGI A,X,NR,NV, IDO,DET,NACT) oooor 

DIMENSION A ( NACT ,NACT) ,X(NACT ,NACT) 00003 

1 F ( IDO) 1,2,1 • • 00002 

1 DO 3 1=1, NR 0000: 

DO 4 J = 1 ,NR 0000' 

4 X(I»J)=0.0 00005 

3 X ( I , I ) =1 ,0 0000( 

NV = NR 0000'. 

2 DET=1 . 0 x 0000? 

Nr i =nr- i . oooo c 

DO 5 K = 1 ,NR 1 * • 0001? 

I R 1 =K+ 1 00013 

P I VOT=0 ,0 00012 

DO 6 I =K»NR 0001: 

Z=ABSF( A ( I , K ) 0001' 

- IF(Z-PIVOT) 6,6,7 - 00015 
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7 P I VOT=Z 
I PR = I 

6 CONTINUE 

IF(PIVOT) 8,9.8 
9 DET=0.0 
RETURN 

8 IF(IPR-K) 10,11,10 

10 DO 12 J = K » NR 
Z = A ( IPR,J ) 

A ( I PR » J ) =A ( K * J ) 

12 A ( K » J ) =Z 

DO 13 J = 1 , NV 
Z = X ( I PR ,J ) 

X( I PR , J ) -X ( K , J ) 

13 X ( K » J ) =Z 
DET = -DET 

11 DET = DET*A ( K. » K ) 

PIV0T=1.0/A(K»K) 

DO 14 J = I R 1 , NR 
A(K,J)=A(K,J)*PIVOT 
DO 14 I=IR1,NR 

14 A( I , J) =A( I , J )-A( I ,K)*A (K, J ) 

DO 5 J= 1 , NV 

I F ( X ( K , J ) ) 15,5,15 

15 X(K,J)=X(K,J)*PIVOT 
DO 16 I = I R 1 , NR 

16 X( I » J ) =X ( I , J ) — A ( I , < ) *X ( K » J ) 

5 CONTINUE 

I F ( A ( N R , N R ) ) 17,9,17 

17 DET=DET*A(NR,NR) 

PIV0T=1.0/A(NR,NR) 

DO 18 J = 1 ,NV 
X(NR»J)=X(NR»J)*PI VOT 
DO 18 K=1 »NR1 
I=NR-K 

SUM=0 • 0 
DO 19 L= I , NR 1 

19 SUM=SUM+A (I,L+1)*X(L+1,J) 

18 X( I , J ) = X ( I , J ) -SUM 
END 

SUBROUTINE MATRA ( M , N , A ,B , MD » ND ) 
Dimension a ( md > nd > , b ( md ,nd j 
DO 3 -1 = 1, M 
DO 3 J = 1 » N 
3 B ( I , J > = A ( J , I ) 

RETURN 

END 

SUBROUTINE MAMUL (M,N,L,A,B,C, MD , ND , LD ) 
DIMENSION A(MD,ND) »B(ND,LD) , C ( MD » LD ) 

DO 2 I = 1 , M 
DO 2 J = 1 > L 
5=0. 

DO 1 K= 1 » N 

1 S=S+A ( I »K ) *8 ( K > J ) 

2 C( I ,J)=S 
RETURN 
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FINIS 
EXECUTE, 3. 

• 4 

, o.nonooonoooE on 2. oooononnooE 

1 .OOOOOOOOOOE 00 0. OOOOOOOOOOE 

' 4. OOOOOOOOOOF 00 3. OOOOOOOOOOE 

5. OOOOOOOOOOE 00 4 . OOOOOOOOOOE 
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0. OOOOOOOOOOE 00 3. ooonopoonoE 
3. OOOOOOOOOOE 00 0 . OOOOOOOOOOE 
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PROGRAM REVEIG2 

A short third program, REVEIG2, was formulated to carry out the cal- 
culation in equation (3-5) of Appendix 3, and to compute the corresponding 
matrix K(U') from equation (8) of the text. The program printout notation 
requires special description however* The input to REVEIG2 was an arbitrary 
real diagonal matrix, G, and an orthonormal matrix X 1 *, as in equation (3~5). 
The latter matrix was, however, the result obtained from REVEIG for a sym- 
metric matrix, and called U 1 or UP in REVEIG instead of X'*o It is there- 
fore called UP in REVEIG2 , also, necessitating the introduction of the 
program description UPP, or U 1 ! , for the calculation in equation (3-5)* 

Thus in REVEIG2 the following is calculated 

UPP = G (UP) (7-1) 

and 

K = (I+UPP^a-UPP) (7-2) 
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-COOP, , YOST* S/1S/2S .3,5000. 

-FTN ,L,E. 

| PROGRAM REVE1G2 

DIMENSION UP(25,25) , G ( 2 5 ) , AUX ( 2 5 , 2 5 ) , UPP ( 25 , 2 5 ) , FK ( 2 5 , 2 5 ) 

DIMENSION AUXK25.25) ,AUX2(25,25) 

C e FIND THE K MATRIX DERIVED FROM UPP=G*UP WHERE G = B, TO THE MINUS ONE-HALF 
C AND K = (I+UPP) INVERSE TIMES ( I -UPP ) 

1 READ 899, N 

I F ( N ) 2,2,3 

2 STOP 

3 PRINT 898, N 

898 FORMAT ( 16H 1UP MATRIX* N= I 3 / ) 

899 FORMAT < 15 ) 

DO 48 1=1, N 

READ 897, (UP ( I , J ) ♦ J= 1 * N ) 

48 PRINT 897, (UP( I ,J) ,J=1 ,N) 

897 FORMAT ( 4E20 • 10 ) 

PRINT 896, N 

896 FORMAT ( 1 5H0G MATRIX, N=I3/) 

READ 895, (G( J ) ,J = 1,N) 

47 PRINT 895, (G( J),J=1,N) 

895 FORMAT (4E20.10) 

C 

DO 46 I = 1,N 
DO 46 J = 1 , N 

46 UPP ( I , J ) = UP ( I , J ) * G ( I ) 

PRINT 894 
DO 45 I = 1 , N 

45 PRINT 893,1 , ( UPP ( I , J ) , J = 1 , N ) 

894 FORMAT (11H0UPP MATRIX/) 

893 FORMAT ( I 5/ ( 1X7E17 . 10 ) ) 

C 

DO 44 I = 1 , N 
DO 44 J = 1 , N 
IF(I-J) 43,42,43 

42 AUX ( I , J ) = UPP( I , J ) + 1. 

AUX2 ( I ,J)=1.-UPP( I ,J) 

* GO TO 44 

43 AUX ( I , J ) = UPP( I , J ) 

AUX2 ( I , J ) =UPP ( I , J ) 

44 CONTINUE ' 

CALL MATALG ( AUX , AUX 1 , N» N , 1 » DE T , 25 ) 

CALL MAMUL ( N , N , N , AUX 1 , AUX2 , FK , 25 , 2 5 , 2 5 ) 

C 

PRINT 892 

892 FORMAT ( 9H0K MATRIX/) 

DO 40 I = 1,N 

40 PRINT 891, I, ( FK ( I , J ) , J=1*N) 

891 FORMAT ( I 5/ ( 1X7E17 . 10 ) ) 

c i 

DO 39 1=1, N 
3 9 G( I ) = 1./G( I )**2 

PRINT 890, ( G ( I ) , 1=1 ,N) 

890 FORMAT ( 9H0B MATRIX / ( 1 X7F 1 7 . 10 ) ) 

GO TO 1 
END 
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SUBROUTINE MATALG ( A . X .NR » NV » I DO , DET , NACT ) 
DIMENSION A (NACT .NACT) , X ( NAC T , N ACT ) 

I F ( I DO ) 1,2.1 

1 DO 3 1 = 1, NR 

DO 4 J = 1 » NR 

4 X ( I * J ) =0 . 0 
3 X ( I , I > = 1 . 0 

NV = NR 

2 DET =1.0 
NR 1 =NR- 1 

DO 5 K = 1 » NR 1 
I R 1 =K+ 1 
PI VOT=0 .0 
DO 6 I =K , NR 
Z=ABSF ( A ( I ,K) ) 

IF(Z-PIVOT) 6,6,7 

7 P I VOT= Z 
I PR = I 

6 CONTINUE 

IF(PIVOT) 8,9,8 
9 DET =0 . 0 
RETURN 

8 IF(IPR-K) 10,11,10 

10 DO 12 J =K , NR 
Z = A( I PR , J ) 

A ( I PR ♦ J ) =A ( K » J ) 

12 A ( K » J ) = Z 

DO 13 J = 1 , NV 
Z = X ( I PR , J ) 

X( IPR,J)=X(K,J) 

13 X ( K , J ) =Z 
DET =-DFT 

11 DET = DET*A ( K » K ) 

P I VOT =1.0/A(K,K) 

DO 14 J= I R 1 , NR 
A(K,J)=A(K, J)*PIVOT 
DO 14 I = I R 1 , NR 

14 A( I ,J)=A( I , J ) —A ( I ,K)*A(K,J) 

. DO 5 J = 1 ,NV 

I F ( X ( K » J ) ) 15,5,15 

15 X<K,J)=X(K, J)*PIVOT . ' 

DO. 16 I = I R 1 , NR 

16 X( I , J ) = X < I , J ) —A ( I > K ) *X ( K » J ) 

5 CONTINUE 

I F ( A ( NR ,NR ) ) 17,9,17 

17 DFT=DET*A(NR,NR ) . . 

PIVOT =1»0/A(NR, NR) 

DO 18 J= 1 , NV 

X ( NR , J ) =X(NR,J)*PIVOT 

DO 18 K=1 ,NR1 

I =NR-K 

SUM=0. 0 

DO 19 L= I ,NR 1 

19 SUM=SUM+A( I ,L+1)*X(L+1 ,J) 

18 X ( I * J ) =X ( I , J ) - SUM 

END 42 
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0000 ! 
0000 
0000 . 
0000 
0000 ' 
OOOO 1 
OOOOl 
0000' 
OOOOl 
0000' 
000 1 < 
0001 
0001 . 
0001 
0001' 
0001 ' 
OOOlf 
0001 ' 
000 II 
0001 * 
0002 ( 
00 02 ' 
0002; 
0002 ' 
0002 < 
000 21 
00021 
0002 ’ 
00021 
0002 ' 
0003C 

0003; 

0003, 

0003' 

000 2> L 

000 3 : 

000 3 ^ 

0003’ 

00031 

0003' 

0004( 

0004-; 

0004; 

0004; 

0004^ 

0004! 

0004( 

0004' 

00041 

0004 

000 51 

0005! 

00051 

0005; 

0005/ 

0005! 



SUBROUTINE MAMUL ( M , N , L , A , B , C ,M k D , ND ♦ LD ) 
DIMENSION A(MD,ND) ,3(ND,LD) ,C(MD,LD) 

DO 2 1=1, M 
DO 2 J=1*L 
S= 0 • 

DO 1 K=1,N 

1 S=S+A( I ,K)#B(K,J) 

2 C( I » J) =S 
RETURN 
END 



END 

FINIS 

EXECUTE. 

4 

9.4067179924E-02 
3. 1644217545E-01 
5. 5556147215E-01 
7.6312984866E-01 
2.oooonooonoE 00 



-2. 3.7 54675249 E-01 
-6.2782339524E-01 
-4. 3075473215E-01 
6. 03207829 4 OE-01 
4.0000oo^>oooe on 



3.7708598062E-01 
5.6007619310E-01 
-7.0046990430E-01 
2.3122009925E-01 
6 . OOOOOOOOOOE 00 



i 



-8.9024101645F- 
4.3819719533E- 
-1 .2305981130E- 
1 .7619122150E- 

8 . nnnnnnnnnoc 



i 
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APPENDIX 8 



NUMERICAL RESULTS 

The following 3x3 unsymmetric matrix was tested in REVEIG so that 
the results could be compared with the known eigenvalues and eigenvectors o 
This example was taken from Lanczos f2]. 



33 


16 


72 


=24 


-10 


-57 


-8 


-4 


-17 



Complete Eigenvalue Analysis 



X = 1 


X - 2 


X = 3 


X = 1 


X = 2 


X * 3 


U 1 


u 2 


u 3 


V 1 


v 2 


v 3 


-15/4 


-4 


-4 


1/4 


0 


4/3 


3 


13/4 


3 


0 


-1/3 


4/3 


1 


1 


1 


1 


1 


1 



The computer results compare exactly with the above answers as 
shown in the accompanying data sheets. 

The following 4x4 symmetric matrix was tested in REVEIG so that 
the results could be compared with its known eigenvalues and eigenvectors. 



2.8 


-0.8 


0.0 


0.0 


-0.8 


1.4 


-0,6 


0.0 


0.0 


-0.6 


1.0 


-0.4 


0.0 


0.0 


.--0.4 


0.4 



Xi » .108 


X2 * -685 


X 3 = 1.612 


X 4 = 3.194 


X 1 


x 2 


x 3 


x 4 


.0940 


-.2375 


.3771 


-.8902 


.3164 


-.6278 


.5601 


.4382 


.5556 


-.4308 


-.7005 


-.1231 


.7631 


.6302 


.2312 


.0176 
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An additional 3x3 symmetric A was tested in order to verify the 
program and the theory for 3x3. These results are also shown, and they 
contain the same UP and VP as in REVEIGl for the same case. 
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52 553 010 1 0 _ 30 
52513 Q3C10C40 
52513 Q 3 G 1 0 3 4 0 
52745 Q7GL0GRC 
52771 Q7CL0LLC 
52736 Q7 C FLOAT 
52756 Q7QLDf'C3 
53010 Q7GLDID2 
47514 SCALE 
47650 Q3QXMCDF 
47455 OSCXTOI 
47364 Q2GL0AUA 



52567 0101C53C 
52513 Q3C10140 
52513 G3Q1044C 
52753 07QL0C0C 
53006 Q7QL0CID 
52733 07QL01C2 
52765 Q7QLDCC3 
5 3 0 i 7 Q70L0R02 
47557 CXSORT 
47456 XTO 1 
47377 L0GF 
47347 OV FLU 



• , TRACE = 6.000C000CE CO 

' 'l 



ITERATE 


• 

REAL PART 
5 . 15470054E 


00 


ALMOST TRIANGULAR FCRM 

TRACE = -6 . OOOOOCCOE CC 

LAGUERRE ITERATIONS 
IMAG. PART 

0 


I TERATE 


3 .000 16635E 


00 


* 0 


EIGENVALUE 


3 .OOOOOCC IE 


00 


0 


I TERATE 


2.3331 11560 


00 


0 


ITERATE 


1 .99999999E 


00 


c 


EIGENVALUE 


1 .999 9999 9E 


00 


0 


ITERATE 


1 .OOOOOOOCE 


00 


0 


EIGENVALUE 


1 . COOOOOOOE 


00 


0 



SU MATRIX 



SUM OF EIGENVALUES = 6.00CC0CCCE CO 



1 



-3. 7500000 00 3E 00- 
2 

3. 0000000007 E 00 
; 3 

1 .OOOOOOOOOCE 00 



3. 999999996 dt 00-4 . CCOOOCO 1 04E CO 
3 . 249999995 C L 00 3 . CCOOOCO 1 76E 00 
1 .OOOOOOOOOCE 00 1 .CCOOOCOCOOE 00 






EIGENVALUE 

EIGENVALUE 

EIGENVALUE 



flV MATRIX 



3 .COCCOOOCE 00 
1 . OOOOOOOCE 00 
2.00000000E 00 

SUM OF EIGENVALUES = 



TRACE = 6.0QCOOOCOE CO. 
ALMOST TRIANGULAR FCRM 
TRACE = 6. OOOCCCCCE CO 

0 

. 0 
0 

6 . OOCCOOOCE 00 



1 



2.4999999999E-01- 
- 2 

0 - 

3 

1 .OOOOOOOOOCE 00 



•1 .45519152266-10 1 . 22 333 3332 1 E 00 
•3.3333333352E-0 1 1 ‘.33333233 18E 00 
1 .OOOOOOOOOCE CO l.OOOOOCOOOOE 00 



8W MATRIX 



£-oicrttoGoKj/\L. 



i 



i- 



6.2499999913E-02 

2 

-2. 9103060456 E- 10 
3 

4. 656612873 GE-10 
CU MATRIX 



8. 4 40 1 1 0 832 2 E- 10-2. 6 1934474 1 IE-09 
-8.3333331742E-02-5.C789737523E-09 
-2. 0954757929 E- 09- 3. 33 33 33232 4E -01 



<e ORTHOGONAL 



1 



u 



-3. 750000000. 5 E 00 

. 2 

3 . 0000000007 E 00' 



3 . 999999 996 6 E 00 4 . CCOOOCO 1 04E 00 



•3 . 24 99999950 F 00- 3 . CCOOOCO 1 76E 00 



1 .OOOGOuOOOGE CO-1 .OOOOGCCOOCt- OC-1 . GCCOOCCOCOE 00 
GW MATRIX 



6 



1 

6 . 24999999 1 3 E-0 2- 
2 

-2.9103830456E-10 

3 

4. 656612873 CE- 10 
U-PRIME MATRIX 
1 

- 1 . 5000000C 1 2E 01 
■ 2 

1 .2000000011 E 01- 
3 

4.00C0000028E 00- 
V-PRIME MATRIX 



1 



8. **40 1108322 E- 10 2 . 6 1 934 474 1 1 E-C 9 
0.333333 17a 2 F-0 2 5 . 8789 7 37 523E- C9 
2.0954757925E-09 3 . 3333332324E-01 



1 . 38564C65S 1 E 01 6 . 9282033532E 00 

■1 . 1258330339E C 1 -5 . 1961 52 53 1 8E CC 0RXV\0 t'JOfc.MM- 

3.464 1016482E 00-1 . 73205C8338E 00 



1 .0000000007E 00- 
2 

c- 

3 

4 . 000000CC28 E 00 



•5.0409513526E-10 2. 30940 1 1 096E CO 

■ 1 . 1 54 7 OC 550 1 E 00 2 . 30940 1 1 C 9 1 E 00 OfcTOkJ0&Rl\L 
3.464 101648 2E 00 1 . 73205C8338E 00 
V-PRIME-1 RANSPOSE TIMES U-PRIME 



i 1 

9.9999999951E-01- 

2 

-3. 7252902 935E-09 
3 

3.25962901 12E-09 
K MATRIX 
1 

-9.2698519866E-02 

2 

1 . 0378667465E 00 
3 

4.6349262993E-02- 
KT MATRIX 
1 

-9.2698525495E-02 

2 

1 . 1984252946E 00 
3 

8.02792774 16E -02- 



ITERATE 

ITERATE 

EIGENVALUE 



REAL PART 
3.41990099E 00 
3. 193 9.9497 E 00 

3 . 1937784CE DO 



LAC-UERRE ITERATIONS 
IMAG. PART 

0 

C 

.0 



-1 .210 719347 CE-08 1 . 79 27 9595 6 1 E- 08 
9.999999995 1C-01 3. 5041 0 1 1 870E-08 - I 

1 .2922 100722C-08 1 . CC000C0008E OC ■ 

1 . 1984252959E 00 P . C 279 2932 37E- C2 
5. 1 04508559 3 E-0 1-8. 9881 8 955 36E-01 
-5. 992 1 264054E-0 1 9 . 5986 0 24 1 95E-0 1 

1 .0370667592E CO 4 . 6349 2 62997E- 02 
5. 1 045087465 E-0 1-5. 9921 2647 1 4E-G1 
-8. 9881 09 6 68 4E -01 9 . 5986 036 1 27E-0 1 

TRACE = 5 . 60000000E CC 

ALMOST TRIANGULAR FCRM 
TRACE = 5. 60GCCUCCE CO 



ITERATE 
ITERATE 
I TERATE 

EIGENVALUE 



ITERATE 

ITERATE 

EIGENVALUE 



2 .4555268 RE 00 
1 .6166923 2 E 00 
1.611 7804CE 00 

1.61 1 7804 CE > 00 



1 .03874635E 00 
6 .8564267 6E-0 1 

6 .85642675E-01 



0 

0 

0 

G 

0 

G 

0 



ITERATE 

EIGENVALUE 



1 .0879845 1 E-0 1 
1 .G879045 1 E-0 1 



0 

0 



GU MATRIX 



Mif'i Ut- C IUU^VALUlD “ 



J) • O J U V W 1. 



1 

1 .2326497265E-01-5.93P058 177CE-01 1 . 6308 5294 856 00-5 . 0526970039.. 01 
2 

4. 1466360672E-01-1 .04GSC77692E 00 2 . 422264 30 546 00 2 .487054642 1 8' C 1 

3 

7.28C0385559E-01-7. 1 4 1 066663 5 C-0 1-3 . G2945C75 5 1 E 00-6.9844462315E 00 

4 .. j 

1 .0000000C0CE 00 1 .OCOOOOOOOCE 00 1 .CC000CC000E 00 1 .0000.0000008 CC_ 

TRACE = 5.6OCC0GC0E 00 j 

ALMOST TRIANGULAR FORM j 

TRACE = 5.600000CCE CO j 

LAGUC-RRE ITERATIONS 



ITERATE 


REAL PART 
3.4 1990Q99E 00 


IMAG. PART 

0 


ITERATE 


3.1 939949 7E 00 


0 


EIGENVALUE 


3. 1937784 8E 00 


0 


ITERATE 


2.45552682E 00 


0 ' 


ITERATE 


1 . 6 1 669232 E 00 


0 


ITERATE 


1 .6 1 1 78040E 00 


0 


EIGENVALUE 


1 .61178040E 00 


0 


ITERATE 


1 .0387463 5E 00 


/N 

u 


ITERATE 


6 . 8 5642675 E- 01 


0 


EIGENVALUE 


6. 856426 7 5E -01 


c 


ITERATE 


1 .0879845 1 E-0 1 


0 


EIGENVALUE 


1 .08798451 E-0 1 


0 



SUM OF E-I GENVALUES = 5.60CC00CCE CO 

9 V MATRIX 

1 —-\ 

1 .2326497265E-01-3.938058 1770E-0 1 1 . 62CB 5294 85E CC-5. C526979089E Cl 
2 

; 4. 1466360672E-01-1 .0408077692E CO 2 . 4222 643C54E 00 2 . 487054642 1 E Cl 

3 

7.2800385559E-C1-7. 1 4 1 066663 5E-0 1 -3 . 0294 5075 51 E OC-6 . 984446251 58 CC 

4 , 

L.1 .OOOOOOCOOOE 7 00 l.OOOOOOOOOCC 00 1 . OOOOOCCCOOE 00 1 .OOCOCOOCCOE OC 
QW MATRIX — * 

( ’ - 

1 

1 .7171297740E 00 2 . 674642 0 1 9C E-08 6 . 72462636 77E-08- 1 . 0477378964E-C9 
2 

2.6746420190E-08 2 . 7485 1 2 1 65 4 E 00 9 . 26665961 74E-08-7. 56699591 88E-C9 

3 

6. 7346263677 E-08 9 . 2 666596 1 7 4 E-OS 1 . 87046 1.7582E 0 1 -6 . 5 1 92580 22 3E-C 8 

4 

-1 . 0477378964 E-09-7. 566995 9 188E-09- 6. 5 1 92 580223E- C8 3 . 22 1 302 1 Q44E 03 
U- PR I ME MATRIX 
1 

9. 40671 7 9924 E-02-2. 375467 524 9L-01 3 . 7703 5980 62E- 01-8. 9024101 645E-C 1 

3. 1 6442 1754^E-0 1-6. 2782 33 952 4 8-01 5.CC076 19310E-01 4 . 38 1 97 1 9533E-C 1 
3 

5. 55 56 1 472 15E-0 1-4. 307 547 321 58-0 1 -7 . 00 V 46-9904 30E-0 1- 1 . 230598 1 1 3 OL-C 1 
1 4 

7. 631 2984866 E— 0 1 6 . 0 320 78 294C E-0 1 2 . 3 1 220C99 25E-0 1 1 . 76 1 9 1 22 1 50E-02 






> 

i 

i 

i 

! 



I 

. 

: 

I 



l 



! 

i 

i 

. 



V-PRIME MATRIX 



1 



9. 406 7 179924 E-0 2-2. 3 754 675 24 9 E-0 1 3 . 77C35 98C62E-0 1- 8 . 9024 10164 5C--0 1 
2 

3. 164 42 17 545E-0 1-6. 2782339 52 4 E-0 1 5 . 6007 6 1 93 1 0E-01 4.38197195 338-01 



V-PRIME-TRANSPCSE TIKES U-PRIKC 
1 



1 .OOOOOOOOOCE 00 1 .22140024 1142-08 
2 

1 . 2340024 1 1 EE-03 9 . 999999999 5 E-0 V 

1. 1 8889 1 4 742 E-06 1 . 29330 1 4650 E-03 
4 

-1 .0913936421 E-l 1 -7 . 058078 640 6 E-l 1- 
K MATRIX 
1 

-4. 001 7766878 E-09-1 .499670 178 dF 00 
2 

1.49967012126 00 1 . 8 7283 1 4 890 E-03 

3 

-3.397041 1444E-01 4 . 0592365 2 1 2E 00 

4 

-1.561 68074 15 E 00-3.904 622 333 2E-0 1 
KT MATRIX 
1 

4.001 7766873E-09 1 . 499670 1 736‘E 00 
2 

-1 . 4 99670 121 2 E 00- 1 . 3 720 3 1 4 89 OE-O 8 

3 

3. 39 7041 1 444 E-01-4. 0 5923652 I 5E 00 

4 

1 .561680741 5E 00 3 . 904622 3332E-0 1 



1. 180891 4 742E-08-1 . 09 1 59 36 4 2 1 t- 1 1 
1 .29330 146 59E-08-7 . 958078 6 4C6E- 1 1 
1 .0000 OOCC COE 00-2.6886937121 E-l C 
2.6886937121E-1C 1 . COOOOOOOCOE £C 

3.59704 11263E-01 1 . 561 680759 1 E CO 
4.0592364591E 00 3 . 90462 1 7 602E-C 1 

j 

3. 2716343412 E-CS-1. 92 42 02029 It CO 
1 .9242020085E 00 1 . 8699 2 1 1 068E-C E 

3.597041 1263E-0 1-1 .56 16307591 E CC 
4. 0 592 3 64 5 91 E CO-3.904621 76C2L-C1 
3.2716343412E-08 1 . 9242020291 1: CC 
1 . 92 42 0 200 85 E CO- 1 .869921 1 0682- C 8 



iun up ciuciivncucj 



-J • CWVVW Wfc 



U MATRIX 
1 

3.07724 17209E-01 
2 

9.00121 16814 E-0 1 
3 

1 .0000000000E 00 



ITERATE 

ITERATE 

EIGENVALUE 



ITERATE 

ITERATE 

EIGENVALUE 



ITERATE 

ITERATE 

EIGENVALUE 



V MATRIX 
1 

3. 07724172 09 E-0 1 
2 

9.00121 1 68 14 E-0 1 
, 3 

1 .0000000000E 00 
W MATRIX 
1 

1 .9049122834E 00 
2 

2. 9103830457 E- IV 
3 

1 .3387762010E-09 
U-PRIME MATRIX 

'• 1 

2. 22958 54 4 60 E-01 
2 

6.5217400456E-01 

x 3 

7. 245 40 23707 E-01 
V-PRIME MATRIX 
1 

2 .22958 54 4 60 E-01 
2 

6. 521 74 00 456 E-01 
3 

7. 24 5 40 23707 E-01 
V-PRIME-TRANSPOSE 
1 

9. 999 9 999998 E-0 1 
2 : 

3.6379788071E-11 

3 

1. 1641532182E-10 



5.8228119506E-0 1 7. 4412236765E 00 
9. 1 189723 152E-01-3. 6548905973E 00 
1.000000000CE 00 1.0000000000E 00 



JRACE = 5 . 2000000CE 00 

ALMOST TRIANGULAR KCRM 

! ( 

TRACE - 5 . 2000000QE 00 



REAL PART 


LAGUERRE ITERATIONS 
IMAG.. PART 

0 


3 . 322 1 8.33 7E 00 


3.1 9296062E 00 


. 1 
0 


3.19293436E 00 


2 . 16568534E 00 


, 0 


1 .547 13834E 00 


* o •/ 


1.54713834E 00 


0 


7 .69321 492E-0 1 


0 


4 • 5992729 9E-0 1 


0 


4 . 59927299E-0 1 


0 



SUM OF EIGENVALUES = 5.20000000E 00 



5.8228 1 19506E-0 1 7 . 44 122367 65E 00 
9. 1 189723152E-01-3.6548905973E 00 
1 .000000000CE 00 1.0000000000E 00 



2. 910 38 30 457 E-l 1 1 . 33877620 10E- 09 
2.1 70607951 CE 00 2 . 0954 757929E-09 
2.0 95 475792 9E -09 6. 9730035083E 01 



3.9 5222072 1 4 E-0 1 8. 91 11 636088E-01 
6. 1894 94 8008 E-0 1-4. 37687798 41 E-01 
6. 787491 601 9E-01 1 . 1 9754008 16E-01 



VP 



-3.952228721 4 E-01 8. 9111636088E-01 
-6. 189494800 8E-0 1-4. 37687798 4 IE- 01 ** 
6.787491 601 5 E-01 1. 1 9754,008 1 6E-0 1 
TIMES U-PRIME 



UP 



3. 637978807 IE- 1 1 1 . 1 641 532 1 82E- 10 
1.000000000CE 00 1.74622982 74E- 1 0 
1. 7462298274 E- 10 1 . 00000CCQ00E 00 



T 



-1.2732925 82 5 E-l 1 1 . 4471 543472E 00-2 . 30 15 283608E-0 1 
1 2 

-1 .4471543472E 00 2 .'7648638932E-1 0 1 . 5425447924E 00 
2. 301 5283620 E- 01-1. 5425447925E 00-1 .6734702512E-10 
KT MATRIX 
1 

1 .2732925825E-11-1 .4471543472E 00 2.3015283608E-01 
2 

1 .4471543472E 00-2.7648 638933E- 1 0- 1 . 5425447924E 00 
3 

-2. 3015283620 E-01 1 . 5425447925E 00 1 . 67347025 1 2E* 1 0 



The K matrix derived from REVEIG using a 4 x 4 symmetric A was used 
as input into REVEIGl. The results for UP and VP from REVEIGl compared 
with those of REVEIG. In addition, UP equalled VP as required in this 
case. 

Other sample results are included, as discussed in Section 4 of the 

text. 
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K MATRIX, N= 4 = 

-4.001 7766878E-09 
1 . 4 99o70 1 2 1 2E 00 
-3.397041 1 444E-0 1 
-1.561 68074 15E 00 



W * 9 YM M ETRIC. 

-1 .49967i. 17G6E 00 
1 .87283 I4899E-08 
4.0592 3652 1 3E 00 
-3.90462233828-01 



3 . 397041 1 263E- C 1 
-4.0592364591E 00 
3.27163434 1 2E-08 
1 .92420200858 00 



1.561680759 

3.9C4621760 

■1.924202029 

1.869921106 



U-PRIME MATRIX - V PR-tKE: 
1 



,3754 675252 8-01 3 . 7708598059E-0 1 -8 . 90241 01 64 5E-<)1 
,27823395268-01 5.60076 1931 4E-01 4 . 38 1 971 9532E-Cj1 
, 307547321 3E-0 1-7. 00469904 33E-01-1 .2305981 130E-01 
,03207829428-01 2 . 5 1 220C9924E-0 1 1 . 761 91 22 1 84E-b2 
V-PRIME MATRIX = TJ PfMMF 



9.40671 79903E-02-2, 
2 

3.1644217541 E-01-6, 

3 

5. 555614721GE-01-4, 

4 

7. 63 12984867 E-01 6. 



1 



9. 4067 1 78376 E-02-2 .37546758628-0 1 
2 

3. 164421 7654 E-01-6. 27823406358-01 

3 

5. 5556 1 48578 E-01 -4. 3075472992E-0 1- 

4 

7. 631 2983 85 2 E-01 6.0 32078 1 70 2E-0 1 
VPT X VP 
1 

1.0000000001 E 00-1 .2267264538E-08 
2 

-1.2267264538E-08 1 .000000000CE 00 

3 

-1.185981091 1 E-08-1 .295848051 IE-08 

4 

8. 1854523157 E- 12 3.81987774748-1 1 
VP X VPT 
1 

1 .0000000018E 00 4.3291947804E-09 
2 

4. 32919478048-09 1 . 0000000099E 00 

3 

5. 329 63 8 9523 E- 10 2. 249 1 803975E-09 

4 

-4. 453 7955546 E-09-4. 844 878276 2E-0 9 
UPT X UP 



3. 7708 59 82 32 E-01-8. 902 4101 633E-C1 
5. 60076 19757E-01 4 . 38 1 97 1 9545E-01 
7. 0046990532E-0 1-1 .2305981 148E-01 
2.31220C8239E-01 1 .76 19 1 22242E-02 

1 . 1 85981 09 11E-08 8 . 18545231 57E-1 2 
1.2958 48051 IE- 08 3 . 8 1 98777474E-1 1 
9.9999999993E-01 2.501 1 104298E-10 
2.5C11 1 04298E- 1 0 9. 9999999998E-0 1 

5.3296389523E-10-4.4537955546E-09 

2.24918C3975E-09-4.8448782763E-C9 

i 

1 . 0000000073E 00 1 . 1 539555089E-08 
1 . 1539555089E-08 9 . 9999998095E-C 1 



1 

T .000000000 1 E 00 
2 

1 .2369127944E-08 

3 

1 . 1797965271 E-08 

4 

5.91 17155614 E- 12- 
UP X UPT 

j 

1 

9.999999981 1 E-0 1 • 
. 2 

-4. 3364707 38 1 E-09 

3 

-5. 4933479986 E-10- 

4 

4.41 1 0493036E-09 
VPT X UP 

"i 1 

1 .0000000001 E 00 
2 

5.820766091 3 E- 11 

3 

-4. 7293724492 E- 11- 

4 

8. 6401996666 E- 12 



1 .2369127944E-08 1 . 1797965271 E-08 5 . 9 1 17 1 5561 4E- 1 2 
l.OOOOOOOOOCE 00 1 .2914824765E-08-4. 0927261 580E-1 1 
1 .2914824765E-08 1 . 0000000001 E OC-2. 1 941 559680E-1 0 
•4.092726 1S8CE-1 1-2. 1941559680 E-10 1 . OOOOOOOOOOE 00 

•4. 336470738 IE-09-5. 49334 79986E- 1.0 4.411 049303 6E-09 
9.999999901 3E-0 1-2 . 3073880584E-09 4.81531 96985E-09 
■2. 307388058 4 E-09 9. 9999999269E-0 1- 1 . 1431 325220E-08 
4.815319698SE-09-1 . 1 431325220E-08 1 . 0000000 1 92E 00 

4.36557456850-1 1- 1 . 09 1 393642 1 E- 1 1 9 . 09494701 75E- 1 2 
l.OOOOOOOOOCE 00-2.5465 8516 50E-1 1 3 . 7744030 1 24E- 1 1 
■3.2741 80 9264E-1 1 1 . OOOOOOOOOOE 00 5. 68434 1 8860E-1 2 

-3.092281986CE-1 1 2 . 50 1 1 1 04299E- 1 1 9 . 9999999998E-0 1 

53 > * " 



I 



K MATRIX, N= 3 - 5 <SVJ- 5 Y H I C 

-1.27 32925825 £-1 1 1 .447 1 543472E 00 -2 . 30 1 5 283608E-0 1 

-1.4471543472E 00 2 . 76486389 33E- 1 0 1 . 542S447924E 00 

2 .301 52836 20E-0 1 - 1 .5425447925E CC - 1 . 67347025 1 2E- 1 0 

U-PRIME MATRIX r V PdlKt' 

1 

2. 2295854460E-01-3. 952228721 EE-01 8 . 9 1 1 1 636085E-0 1 
2 

6. 52 1 7400455 E-p 1-6. 189 494 8013 E-0 1-4. 3768779843E-01 
3 r 

7. 2454023702 E-Q1 6. 787491 60 1 9E-0 1 1 . 19754008 12E-01 
V-PRIME MATRIX * V PM HF 



1 

2.2295854451 E- 
2 

6.521 7400463E- 
3 

7. 2454023706 E- 

VPT X VP 
1 

1 .000000000CE 
2 

-1. 4551915228 E- 
3 

-1 .2732925825E- 

VP X VPT 
1 

1 .OOOOOOOOOOE 
2 

-2. 1 82 7 87 2843 E- 
3 

-1 . 23691 27944 E- 

UPT X UP 
1 

9.9999999991 E- 
2 

-2.91 03830457 E- 
3 

7.821 6544352E- 

UP X UPT 
1 

9.9999999988 E- 

I 2 

7 .2759576 141 E- 
3 

1.091 3936421 E- 

VPT X UP 
1 

9. 9999999995 E- 
2 

-7.2759576141 E- 
3 

-7.275957614 lE- 



01 -3. 9 522287 22 5E-0 1 8.91 11 636090E-01 
01-6. 1 894947 997E-0 1-4.37 68 7798 39E-01 
01 6.787491 601 2E-01 1 . 1975400796E-01 

00-1 .455191 5226E-1 1-1 . 2732925825E- 1 0 
•11 9 .9 99 9 99998 8 E-0 1-1 . 4006 2 1 8407E- 1 0 

•10-1 .4006218407E-10 1 . OOOOOOOOOOE 00 

\ 

00-2. 1827872 843 E-l 1- 1 . 2369 1 27944E- 1 0 
■11 9. 9 99 9 999 99 8 E-0 1 1 . 3733369996E- 1 0 
•10 1.3733369996E-10 9. 9999999991E-01 

-01-2.9 10 3830457E-1 1 7. 82 16 544352E- 1 1 
•1 1 1 . 000000000 1 E 00 1 .7462298274E-10 
•11 1 .7462298274E-1 0 9. 9999999993E-01 

•01 7.2759576141E-12 1 . 09 13936421 E- 1 0 u 
•12 1 .000000000 IE 00-1 . 69 1660 1 4 53E- 1 0 
•10-1. 6916601453 E- 10 1 . OOOOOOOOOOE 00 



-01-3.6379788C71E-1 1-3. 4560798667E- 1 1 
-12 9. 9999 99999 EE-0 1 5 . 4569 6e2 1 05E- 1 2 
■12 2.91C3830457E-1 1 9 . 9999999998E-0 1 



= T 



K MATRIX, N= 4 - 

0 -1 .OOOOCOCCOOE 00 

1 .OOOOCOOC OOE 00 0 

4 . OGOGGOOC OCE O'O 3 . 0000 COOOOOE CO 

5 . OOOOOOOGOOE 00 4 .OOOOCCOOOOE 00 

U-PRIME MATRIX S V 



-4. OOOOOCOCOCE 00 
-3.0000000C0CE 00 

0 

3 . COOOOOOOOOE 00 



-5. OOOOCOOCC 
-4 . CCCOCOOOC 
-3.000000000 



-1 .35802469 15 E-01-6. 1 728 39505 5 E-0 1 -3 . 9506 1 728 36E-0 1 6 . 6666666665E-0 1 
2 

-9.6296296 294 E-0 1 2 . 5 925 9259 34 C-0 1- 7 . 4C74 07 4 1 02E-02 0 

3 

-1 .975308641 7 E-0 1-7. 1 604938272E-0 1 6. 1 728394998E-02-6 .6666666665E-C 1 

4 

1 .2345679009 E-0 1 1 . 9753086 42 2E-0 1-9 . 1 358024699E-0 1-3 . 3333333330E-C 1 
V-PRIME MATRIX ■» U TfUHE 

-1 .358024 69 15 E-0 1-6. 17283950 5 5E-0 1-3. 9506 172836E-01 6 .6666666665E-0 1 
2 

-9. 6296296294 E-0 1 2 . 5925 925 934 E-0 1 -7 . 4074 074 1 02E-02 . 0 

3 

-1 . 975308641 7E-0 1-7. 1 604938272E-0 1 *6 . 1 728 394998E-02-6 . 6666666 665E-0 1 

4 

1.2345679009E-0 1 1 . 9753086422E-0 1 -9 . 1 358024699E -01 -3 . 3333333330E-0 1 
VPT X VP 

9. 9999999993 E-01 -9. 4 58744898 4 E-l 1 5 . 45696 82 1 06E- 1 1 -1 . 7280399334E-1 1 
2 

-9. 4587448984 E-l 1 9.9999999998E-0 1-4. 7293724492E- 11 4. 1 836756282E- 1 1 

3 

5. 456 96 82106 E- 11-4. 7293724 49 2 E- 1 1 1 .0000000001 E 00 6 .54856 1 8 528E- 1 1 

4 

-1 .7280399334E-1 1 4 . 1 836756282E- 1 1 6 . 54836 1 85 28E- 1 1 9.9999999995E-01 
VP X VPT 

9. 9999999988E-0 1-1 .637090463 2E-1 1-7 . 2759 576 14 1 E- 1 2 3. 637978807 1 E-l 1 
2 

-1.6370904632E-1 1 1 . OOOOOOOOOCE OO-'fc. 60609361 55E- 1 1 8 . 5492 501 967E- 1 1 

3 

-7. 2759576 14 IE- 12-8. 6060936 155E-1 1 9. 9999999998E-0 1 1 . 455 1 915228E-1 1 

4 

3. 637 9788 07 IE- 11 8 . 549250 1 967E-1 1 1 . 4551 9 1 5228E- 1 1 1 . 000000000 IE CC 
UPT X UP 

9. 9999999993E-0 1 -9.458744898 4E-1 1 5 . 4569682 1 06E- 1 1- 1 . 7280399334E- 1 1 
2 

-9.4587448984 E-l 1 9. 9999999 99 8 E-0 1-4. 7293724492E-1 1 4. 1836756282E-1 1 

3 

5. 456 96 82106 E- 11-4. 729372 4 492 E- 11 1 . 0000000001 E 00 6 . 548361 8528E- 1 1 

4 

-1. 7 28 0399 334 E- 11 4. 1 836756282E-11 6. 54836 18528E- 1 1 9 . 9999999995E-0 1 
UP X UPT 
1 

9.9999999988E-01-1.6370904632E-1 1-7 . 2759576 1 4 1 E- 1 2 3 . 63797 88 071 E-l 1 
. 2 

-1 . 6370904632E-11 1 . OOOOOOOOOOE 00-8 . 6060936 1 55E- 1 1 8.5492501 967E-1 1 

3 

-7.2759576141 E-l 2-8 . 60609361 55E- 1 1 9. 9999999998E-01 1 .455 1915228E-1 1 

4 • ' 

3. 6379788071 E-l 1 8. 549250 1 96 7E-1 1 1 . 4551 9..1 5228E- 1 1 1 . 000000000 1 E OC 



:-iTl 



VPT X UP 

9. 9999999993 E-0 1-9. 4 58744898 4 E-l 1 5 . 4569682 1 06E- 11-1. 7280399334E- 
2 

?9. 458744 8984 E-l 1 9. 999999999 8E-0 1-4. 72937 24492E-1 1 4 . 1 836756282E- 1 1 
5. 4 569662 106 E-l 1-4. 729372 449 2E-1 V 1 . 0000000001 E 00 6. 548361 8528E-1 1 

-1 .7280399334E-1 1 4. 1 836756282E-1 1 6 . 54 836 1 8528E- 1 1 9 . 9999999995E-0 1 

— 

- 55 



=1 



K MATRIX, N= 4 » 
9 . 0000000000 E 00 
1 .OOOOOOOOOOE 00 
4 . 0 000000000 E 00 
5. OOOOOOOOOOE 00 



uW&rfAMETCUC. 

2. OOOOOOOOOOE 00 
6 . OOOOOOOOOOE 00 
3. OOOOOOOOOOE 00 
4. OOOOOOOOOOE CO 



t 



5. OOOOOOOOOOE 00 
2. OOOOOOOOOOE 00 
7. OOOOOOOOOOE 00 
3. OOOOOOOCCOE 00 



1 . 000000000 

2. OOOOCOOCO! 

3 . OOOOOOOCO. 
4 . OCOOOOOOOj 



U-PRIME MATRIX 



-7. 63440860 18 E-0 1-3. 4408602 157E-0 2-1 . 63440860 19E^-01 6.451 61 29030E-C2 
2 ' 
5.017921 1473E-02-6.315412186 2 E-0 1-8 . 31 54 1 2 1875E-02- 1 . 07526881 72E-j01 

3 1 

-4^301 0752691 E-02-3 . 0 107526872E-02-6. 4301 0752 65E-01- 1 . 935483871 2E-jp 1 

4 

-2. 50 89605734 E-0 1-2. 422 939068 2 E-0 1 1 . 5770609340E-02-4 .6236559 1 40E-0 1 



V-PRIME MATRIX 



-1. 206896551 7E 00-2.0689655 17 5E-0 1-1 .3793 103452E-01 7. 586206897 3E-01 
2 

1 . 7733990 149 E-0 1-1 .6798029557E 00-1 . 6748768472E-01 7.7832512329E-0 1 

3 • 

2.7586206894 E-0 1 2. 7586206902E-0 1- 1 . 48275 86206 E 00-3.4482758626E-0 1 

4 

-3. 2512315277 E-01 2. 463054 1 874E-0 1 6. 4039408874E-01-2.0935960592E 00 
VPT X VP 
1 

1 . 6698536 727 E 00-5 . 2 17 306 89 9 7 E-0 2-4 . 8047756554E- 01- 1 . 9 1 99 689376 E-0 1 

} f > 

1 2 

-5.21 73068997E-02 3.001 3 103936E 00 5 . 8579436492E-02-2 . 075 1 777529E CO 
' 3 

-4. 8047756554 E-0 1 5.8579436492E-02 2. 6557548108E 00-1 .0644276738E 00 
4 

- 1.9199689 37 6E -01-2. 0751 77 7529F 00- 1 . 0644276738E 00 5 .6833458720E OC 



VP X VPT 
1 

2 . 0939357908 E 00 7 .470698 1 50 7E-0 1 -4 . 47086801 44E- 01- 1 . 335 1 452355E OC 
, 2 

7 . 4706981 507 E-01 3. 4870295327E 00-4. 3451 673200E-0 1 -2 . 208 1 5841 2 5E 00 

3 

-4. 4708680 144 E-0 1-4. 3451 67 320CE-01 2 .4696789536E 00-2. 49363003 12E-C1 

4 

-1 .3351 452355E 00-2. 2081 584 1 25E 00-2 . 49363003 1 2E-0 1 4 .9596204719E 00 
UPT X CJP 



1 1 

6.5015865669E-01 5. 6664225799E-02 1 . 443044 1 538E-0 1 6 .9680502568E-02 
2 

5. 6664225799E-02 4 . 5964 10 631 9E--0 1 7 . 3677 367961 E-02 1 .8354337690E-01 

3 

1 .4430441538E-01 7. 367736796 IE-02 4. 4733906290E-01 1 . 1555863876E-0 1 

4 

6.968050256QE-02 1 .8354337690E-0 1 1 . 1 555863876E-0 1 2 . 6696727946E-0 1 
UP X UPT 






j 6. 14901 i 4 459E-0 1-9. 9248 46800 5003 V. 26479 36 1 75E- C 1 1 .674 73696374-01 
2 

-9. 9248468005 E-03 4 . 1 983890238E-0 1 9 . 1 1 36547578E-02 1 .88834161 96E-01 

3 

1 .2647936175E-01 9. 1 1 36547578E-02 4 . 53680 1 9420E-0 1 9.7435541679E-02 

4 

1 .6747369637E-01 1 . 38834 16 1 96E-0 1 9. 743554 1679E-02 3 . 3568582 11 1 E-0 1 
VPT X UP 

1 ~ 

9. 9999999995 E-0 1 2 . 364686224 6E-1 1-3 . 0468072509E- 1 1 2 . 9 1 03830457E-il 1 
2 

-1. 8189894035 E- 12 1 .000000000CE 00-3. 2969 1 82939E- 1 2-2. 7284 84 1 05 3E- VI 
3 

1 .09 1 393642 IE- 1 1-5. 456968 2106 E-l 1 9 . 9999999993E-0 1-2 . 9 1 03830457E-?. 1 1 
K " 

-1 . 455191 5228E-1 1-2.9 103830457E-1 1-3.6379780C71E-11 1 .0000000001 E 00 



K MATRIX, 



N = 



: unsymmetric oj rr+v 'e.E‘R.0 pifl&owh 1 * 



o 



•1.00000000006 00 
4.0000000C00E 00 
5. 0000000 000E 00 

U-PRIME MATRIX 



1 



2.0000000000E CO 
0 

3.0000C00000E 00 
4.0000C00000E 00 



5. OOOOOOOCCOE 00 
2 . 0000000C00E 00 
0 

3. OOOOOOOOOOE 00 



t 



1 . occoooooo 

2. CCC000000 

3.000000000 



-3. 499999999SE 00 4. 2499999997E 00-2.74999999986 00 2 . 2499999999E QO 
2 

3.4999999998E 00-7. 750000000CE 00 4. 2499999999E 00-2. 7499999999E 00 

3 * v 

-4.9999999964 E-0 1 1 . 749999999 5 E 00-2.2499999998E 00 7.4999999986E-01 

4 r 

-1 .2914824766E-10 5.000000001CE-0 1 4. 9999999990E-01-1 . 4999999997E 0.0 



V-PRIME MATRIX 



/ 



-6. 6666666666E-0 1-3. 3333 333337E-0 1-1. 66666 66667E-0 1-1 . 6666666668E-01 
2 ■ 
-3.7500000001 E-01-3.7500000002E-01-2. 8124999999E-01-2. 1 875000001 E-0 1 
3 1 

4. 1666666671 E-0 2-2. 91 66 66 666 9E-0 1-8. 0208 333333E-01 -3 .6458333334E-?01 

*>■2. 9 1666 66667 E-01 4.1 66666 665 6E-0 2-1 . 3541 6666 68E-0 1-6 .979 1666666E-01 

VPT X VP 

6.71 87500002E-01 3 . 3854 166670E-0 1 2.2265625001E-01 3.81 51041 669E-01 
3. 38541 66670E-01 3. 3854 1 66673E-0 1 3. 8932291670E-01 2. 1484375004E-0 1 
2.2265625001 E-01 3 . 893229 1 67CE-0 1 7.6855468751E-01 4.7623697920E-01 
3. 81 5104 1669E-01 2. 1 484375004E-0 1 4. 76236 97920E-01 6. 9563802086E-0 1 
VP X VPT 

6. 1 1 1 1 11 1 1 15E-01 4.5833333338E-01 2 . 6388888892E-01 3 . 1944444447E-C 1 
4.5833333338E-01 4.0820312502E-0 1 3. 9908854 1 68E-0 1 2. 8450520835E-01 
2. 6388888892E-0 1 3.9908854168E-0 1 Q. 630642361 5E-01 3 . 3875868057E-C 1 
3. ^ 944444 447E-01 2.8450520835E-0 1 3. 3875868057E-01 5.9223090279E-01 
UPT X UP 

2.4749999996E 0 1-4. 2874999995E 01 2. 5624999997E 01- 1 . 7874999997E 01 
-4. 2 874 99999 5 E 01 8. 1437499993E 0 1-4 . 83 12499996E 01 3. 1437499997E 01 
2.5624999997E 01-4.8312499996E 01 3.0937499998E 01-2 .0312499997E 0.1 
-1 .7874999997E 01 3. 1 437499997E 0 1-2 . 03 12499997E 01 1 . 5437499997E 01 
UP X UPT , . 









4. 2937499994E 01 -6 . 3062499994E 01 1 . 7062499994E 0 1-2 . 6249999979E 00 
-6,. 30 6 2 49 9 9 94 E 01 9;7937499997E .0 1-2. 6937499993E 01 2.3749999970E 00 
1 .7062499994E 0 1 -2. 6937499993E 01 8 . 9374999970E 00- 1 . 3749999992E CO 
-2.6249999979E 00 2. 374999997CE 00-1 . 3749999992E 00 2.7499999990E 00 

VPT X UP ' 

* ( % % 

9. 9999999991 E-01 ' 2. 582964953CE-1 0-5. 82076609 13E- 1 1-2 .91 03830457E-1 1 
-7.4502774320E-11 1 .0000000003E 00-1 . 9099388737E- 1 1 2. 3646862246E-1 1 
-2. 22617 776 06 E- 10 3. 637978807CE-1 0 9. 9999999993E-01 6 .54836 18528E-11 
-2.2642628509E-1 1 1 .7462298274E-1 0 3. 6379788071E- 1 1 9 . 9999999979E-02 




Only one test is shown for REVEIG2 since the results present no 
particular significance except that the matrix K is not skew- symmetric 
as was anticipated. 
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UP MATRIX, 



4 



9.40671 799 2E-02 
3. 1644217544E-01 
5.55561 47215E-01 
7. 63 12 984 666 E-0 1 

G MATRIX, N= 4 

2. 0C00C00000E 00 

UPP MATRIX 

1 

1 .8813455985E-01-4 
2 

1 .2657687018E 00-2 
3 

3. 33336S8329E 00-2 
'4 

6. 1050387893E 00 4 
!< MATRIX 
1 

9.7463447412 -01 2 
2 

2.874621 9294 E-02 1 

3 

-8.4402690113 -02 8 

4 

1 .5098146321 E-01 7 
B MATRIX 

2 . 4 999999999 E—0 1 6 



I 



-2.3754675249E-C i 
-6.27822 X 9524E-C1 
-4.30754732 14E-01 
6. 032078 29 39 E-01 



3.7708598C62E-01 
5.60076193C9E-01 
-7.C046990430E-01 
2. 3122009924E-01 



-8.9C24 1C164 
4.361971953 
-1 .2305981 13 
1.761912215 



4.OO0OCCQ0OOE 00 



6 . OOOOOCOCOOE 00 8.COOOCOOOOC 



7509350497E-01 7 . 54 17 1 96 1 24E-0 1- 1 . 7804820329E CC 
5 1 1 29358 1 CE 00 2 . 2403047724E 00 1 .752788781 4E $C 
i 5845283929E 00-4 . 2028 1 94258E 00-7 . 383586678CE-G 1 
,825662635 2E 00 1 . 84976C7940E 00 1 . 4095297720E-C 1 

,9933249680 E-0 1 8 . 78359027 1 4E-0 1-2 . 2909850 1 24E-C2 
,0 79978421 4 E-0 1 -7. 1 805 0 1 5088E-0 1 -2 . 772042880 1.E-C2 
,476466871 6 E-0 1 - 1 . 9359 6 353 58E- 0 1 3 . 1 397766744E-04 
,9680847596E-01 2 . 72 15 8096 27 E-0 1 9 . 9224 176882E-0 1 

.250000000CE-02 2.7777777777E-C2 1 . 5625000000E-C2 
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